Maps close to identity and universal maps in the Newhouse domain 



Dmitry Turaev 

Imperial College, London, SW7 2AZ, UK 
dturaev@imperial.ac.uk 

September 7, 2010 



Abstract 

Given an n-dimensional C r -diffeomorphism g, its renormalized iteration is an iteration of g, restricted 
to a certain n-dimensional ball and taken in some C r -coordinates in which the ball acquires radius 1. We 
show that for any r > 1 the renormalized iterations of C r -close to identity maps of an n-dimensional unit 
ball B n (n > 2) form a residual set among all orientation- preserving C -diffeomorphisms B n — > R n . In 
other words, any generic n-dimensional dynamical phenomenon can be obtained by iterations of C r -close 
to identity maps, with the same dimension of the phase space. As an application, we show that any C r - 
generic two-dimensional map which belongs to the Newhouse domain (i.e., it has a wild hyperbolic set, so 
it is not uniformly-hyperbolic, nor uniformly partially-hyperbolic) and which neither contracts, nor expands 
areas, is C r -universal in the sense that its iterations, after an appropriate coordinate transformation, C r - 
approximate every orientation-preserving two-dimensional diffeomorphism arbitrarily well. In particular, 
every such universal map has an infinite set of coexisting hyperbolic attractors and repellers. 

1 Ruelle-Takens problem and universal maps 

A long standing open problem in the theory of dynamical systems is to describe which kind of dynamical 
phenomena can be expected in close to identity maps. It started with a much celebrated paper pQ where it was 
shown that any n-dimensional dynamics can be implemented by a C"-small perturbation of the identity map 
of an n-dimensional torus. The paper seized a lot of attention by physicists, because it proposed a new view 
on the onset of hydrodynamical turbulence; at the same time it caused a lot of criticism. One of the reasons 
for the critique was that the C n -small perturbations constructed in [I] were not small in C n+1 , which is quite 
unphysical. The controversy was resolved in [2] where it was shown that for any r, given any C-difFeomorphism 
F of a closed n-dimensional ball, one can find a C r -close to identity map g of the (n + l)-dimensional closed 
unit ball B n+1 such that the diffeomorphism F coincides with some iteration of the map g restricted to some n- 
dimensional invariant manifold. Thus, the restriction on smoothness of perturbations was removed by sacrificing 
one dimension of the phase space; anyway, other scenarios of the transitions to turbulence had already been 
known. 

From the purely mathematical point of view, the question still remained unsolved: can an arbitrary n- 
dimensional dynamics be obtained by iterations of a C""-close to identity map of B n , i.e. in the same dimension 
of the phase space? The difficulty is that the straightforward construction proposed in [1] does not work for 
high r in principle. Indeed, given an orientation-preserving diffeomorphism F : B n — > R n , one can imbed it into 
a continuous family Ft of the diffeomorphisms such that T\—F and Fq = id. Then, given any N, the map F 
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Figure 1: An illustration to Ruelle-Takens construction. 



can be represented as a superposition of N maps 

F = F N o...oF 1 , where F s = T /N o T^_ i)/N , (1.1) 

that are 0(l/./V)-close to identity. One can then choose N pairwise disjoint small balls D s E B n of radius 
p ~ N^ 1 ^ 71 ^ 1 ^ and define a map cp : B n — > B n such that (j)(x)\ xe D s = x s +i + pF s ( x ~ Xa ) where x s is the center 
of D s ; the positions of the centers are chosen in such a way (Fig.l) that the distances ||x s +i — %s\\ ar e uniformly 
close to zero for all s, hence the map 4> is C°-close to identity. By construction, 4> n \d is linearly conjugate 
to F, i.e. the dynamics of 4> n \d coincides with the dynamics of F. However, the derivatives of </> of order 
k behave as N~ 1 p 1 ~ k ~ TV"- 1 , i.e. at k > n they do not, in general, tend to zero as N — > +00. Thus, an 
arbitrary n-dimensional dynamics can be implemented by iterations of C Tl_1 -close to identity maps of B n , but 
the construction gives no clue of whether the same can be said about the C™-close to identity maps. 

One could try to position the regions D s differently, or make their radii vary, or change their shape. This, 
however, hardly can lead to an essential increase in the maximal order of the derivatives (of (f>) which tend to 
zero as N — > +00. The reason lies in a well-known fact from the averaging theory that the 0((5)-close to identity 
map 

x — x + Sf(x) 

approximates a time shift of a certain autonomous flow with the accuracy 0(5 m ) for an arbitrarily large m (if 
/ G C°°). Hence, the number of iterations necessary in order to obtain a dynamics which is far from that of 
an autonomous flow, has to grow faster than 0(S~ m ), for every m. As we see, in order to obtain such kind of 
dynamics, one has to control a very large number of iterations of close to identity maps, hence decompositions 
much longer than that given by (|1.1|) have to be considered. 

In this paper we propose such a decomposition (Theorem [4J, using which we show that 
an arbitrary C r -generic orientation-preserving n- dimensional dynamics can be obtained by iterations of C r -close 
to identity maps of B n , n>2. 

To make the formulations precise, we borrow some definitions from [3] . Let g be a C r -diffeomorphism of a 
closed n-dimensional ball D. Take an integer m > and any C'-diffeomorphism ip of R n such that ijj(B n ) C D. 
The map g m ,4> = g m is a C r -diffeomorphism that maps B n into R n . We will call the maps g m ,ip 
obtained by this procedure renormalized iterations of g. 
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Theorem 1. In the space of C r -smooth orientation-preserving diffeomorphisms of B n into R n (n > 2) there 
is a residual set S r such that for every map F£c5 r , for every 8 > and for every n-dimensional ball D there 
exists a map g : R n — )• R n , equal to identity outside D, such that \\g — id\\cr < 8 and F is a renormalized 
iteration of g. 

In other words, a generic C r -diffeomorphism B n —¥ R n is, up to a smooth coordinate transformation, an 
iteration of an arbitrarily close to identity map. This theorem is proven in Section 01 Actually, we prove there 
that 

for any 8 > and e > 0, for every orientation-preserving C r -diffeomorphism F : B n — > R n there exists a S-close 
to identity map g, equal to identity outside a given ball D, such that \\F — g m ,ij>\\c r < £ for some renormalized 
iteration g m ,if> = ^ 1 ° g m ° V'l-B™ (it is enough to prove this for one particular ball D that may depend on F, but 
the choice of D should be fixed in advance, not depending on S nor on the accuracy e of the approximation; then 
for other balls the claim will remain true because there always exists an affine conjugacy that takes one ball to 
the other). Moreover, we construct the map g and the coordinate transformation ip in such a way that the balls 
g l (tp(B n )) (i = 0, . . . , m — 1) do not intersect each other — hence, by adding small, localized in g m ~ 1 (ip(B n )), 
perturbations to <?, we may make g m ,%p run an open subset in the space of C r -diffeomorphisms B n — > R n . Thus, 
the set S r (5) of all renormalized iterations of the maps g : D D such that ||<? — jg?||c < 8 contains a subset 
which is open and dense in the space of C-smooth orientation- preserving diffeomorphisms of B n into R n , for 
every 8 > 0. Hence, the intersection S r of these sets over all <5 > is residual (and independent of the choice of 
D), which gives us the theorem. 

The first step in our construction of the approximations of the given map F by renormalized iterations 
is a representation of F as a superposition of a pair of certain special maps and some volume-preserving 
diffeomorphisms (Lemma [1]). Each of the special maps can be realized as a flow map through a kind of 
saddle- node bifurcation (see Fig. 2), reminiscent of the so-called "Iljashenko lips" (see [4]). For the volume- 
preserving diffeomorphisms one may adjust the results obtained in [3] for symplectic diffeomorphisms and prove 
(Lemma [2]) the existence of an arbitrarily good, in the C r -norm on any compact, polynomial approximation 
by a superposition of volume-preserving Henon-like maps. It is known that Henon-like maps often appear as 
rescaled first-return maps near a homoclinic tangency (cf. [SJ[SJ[7])- In this paper we find a kind of homoclinic 
tangency which does incorporate all the Henon-like maps that appear in our volume-preserving polynomial 
approximations. 

Thus, we show that the map F can be approximated arbitrarily well by a superposition of maps related to 
certain homoclinic bifurcations. The last step is to build a close to identity map which displays these bifurcations 
simultaneously. This is achieved by an arbitrarily small perturbation of the time-i5 map of a certain C°° flow 
(the time 8 map of a flow is, obviously, 0(i5)-close to identity). 

Note that the approximation (that we construct in Lemma [2]) of any volume-preserving diffeomorphism 
of a unit ball into R n by a polynomial volume-preserving diffeomorphism is not straightforward, because the 
Jacobian of the approximating diffeomorphism should be equal to 1 everywhere, and this constrain is quite 
strong for polynomial maps. Had the approximation result been true for all volume-preserving maps, i.e. not 
necessarily diffeomorphisms, it would produce a counterexample to the famous "Jacobian conjecture" ; however, 
our approximation uses in an essential way the injectivity of the map that has to be approximated (we represent 
the map as a shift by the orbits of some smooth non-autonomous flow) . 

It should be mentioned that Theorem [1] does not hold true at n — 1. Indeed, if a map F on the interval B 1 
has two fixed points (with the multipliers different from 1), then every close map F has a pair of fixed points 
Pi 2 as well. If such F is a renormalized iteration of a diffeomorphism g, i.e. if F = ip o g m o ip, then g will 
also have a pair of fixed points, ip(Pi) and ip(P2) (at n > 1 this is not true). The interval between Pi and P2 
will therefore be invariant with respect to o g o ip t hence ip -1 0301/1 will be a root of degree m > 1 of the 
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map F on this interval. Now note that the maps of the interval that have a root are not dense in C 2 , according 
to [8]. Thus we obtain that renormalized iterations are not dense either. 

One can check through the proof of Theorem [1] that it holds true for finite-parameter families of orientation- 
preserving diffeomorphisms: 

in the space of k-parameter families F e , e G B k , of C r -smooth orientation-preserving diffeomorphisms of B n 
into R n (n > 2) there is a residual set Skr such that for every F £ € Skr cind for every S > there exists 
g e : D x B k — T D such that \\g E — id\\cr < 8 and F e = ipj 1 o g" 1 o ip s \B-» for some m > and some family ip e of 
C r -diffeomorphisms of R n . 

Thus, any dynamical phenomenon which occurs generically in finite-parameter families of dynamical systems 
can be encountered in maps arbitrarily close to identity (with the same dimension of the phase space) . 

To put the result into a general perspective, we recall that one of the main sources of difficulties in the theory 
of dynamical systems is that structurally stable systems are not dense in the space of all systems [9] [TUJ [11] , 
moreover most natural examples of chaotic dynamics are indeed structurally unstable (like e.g. the famous 
Lorenz attractor [Hj). Understanding the dynamics of systems from the open regions of structural instability 
(i.e. the regions where arbitrarily close to every system there is a system which is not topologically conjugate 
to it) has been the subject of active research for the past four decades. It often happens, and helps a lot, that 
structurally unstable systems may possess certain robust properties, i.e. dynamical properties which are not 
destroyed by small perturbations. For example, systems with Lorenz attractor are pseudohyperbolic (or volume- 
hyperbolic) [12] [T31 rjj] , and this is, in fact, the very property which allowed for a very detailed description of 
them [T51 [TH] [T7] . Another robust property is uniform partial hyperbolicity, a rich theory of systems possessing 
it is actively developing [181 119] . In fact, not so much of robust properties are known, it could even happen 
that beyond the mentioned partial hyperbolicity and volume-hyperbolicity no other robust dynamical properties 
exist. This claim can be demonstrated for various examples of homoclinic bifurcations (see [20] ) . and can be 
used as a guiding principle in the study of bifurcations of systems with a non-trivial dynamics: 
given an n-dimensional system with a compact invariant set that is neither partially- nor volume-hyperbolic, 
every dynamics that is possible in B n should be expected to occur at the bifurcations of this particular system. 
The last statement is not a theorem and it might be not true in some situations, still it gives a useful view 
on global bifurcations. In particular, it was explicitly applied in [21] to Galerkin approximations of damped 
nonlinear wave equations in order to obtain estimates from below on the dimension of attractors in the situation 
where classical methods [55] do not work. 

Theorem [1] gives one more example to the above stated principle: the identity map has no kind of hyperbolic 
structure, neither it contracts nor expands volumes, so it should not be surprising that its bifurcations provide 
an ultimately rich dynamics. 

The same idea can be expressed in somewhat different terms. Let us call the set of all renormalized iterations 
of a map g : D — > D its dynamical conjugacy class. The map will be called C r -universal (cf. [3]) if its dynamical 
conjugacy class is C r -dense among all orientation-preserving C r -diffeomorphisms of the closed unit ball B n into 
R n . By the definition, the dynamics of any single universal map is ultimately complicated and rich, and the 
detailed understanding of it is not simpler than understanding of all diffeomorphisms B n — > R n altogether. 

At the first glance, the mere existence of C-univesal maps of a closed ball is not obvious for sufficiently 
large r. However, Theorem [1] immediately implies the following 

Theorem 2 . For every r > 1, C r -universal diffeomorphisms of a given closed ball D exist arbitrarily close, in 
the C r -metric, to the identity map. 

Proof. Take an arbitrary sequence of pairwise disjoint closed balls Dj C D, a sequence of maps Fj which is 
C r -dense in space of orientation-preserving C r -diffeomorphisms B n — > R n , and a sequence £j —> +0 as j —5- +oo. 
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By Theorem [TJ given any 6, there exist maps gj such that gj is identity outside Dj, some renormalized iteration 
of gj is enclose to Fj, and \\gj — id\\c r < 8. By construction, the map g(x) = gj(x) at x 6 IX, (j = 1,2,...) and 
g(x) = x at x G D\ U^. 1 £>j is C r -universal and <5-close to identity. □ 

As an immediate consequence of Theorem [3J we note that 
arbitrarily close to identity map in the space of C r -diffeomorphisms B n — > B n (any r > 1, any n > 2) there 
exist maps with infinitely many coexisting uniformly-hyperbolic attractors of all possible topological types. 
This is true because a hyperbolic attractor is a structurally stable object: given a map with a uniformly- 
hyperbolic attractor, any C 1 -close map has a hyperbolic attractor topologically conjugate to the original one. 
For every n > 2 there exists a C""-diffeomorphism B n — > R n with a hyperbolic attractor. Hence, infinitely many 
of the maps Fj in the proof of Theorem [5] have a hyperbolic attractor as well. It follows that each universal 
map constructed in Theorem [5] has infinitely many hyperbolic attractors (and repellers, for that matter). 

Thus, we can further pursue the approach of |T| and claim that hyperbolic attractors can be born at the 
third Andronov-Hopf bifurcation. Recall that at the primary Andronov-Hopf bifurcation a limit cycle is born 
from an equilibrium state, and at the secondary Andronov-Hopf bifurcation a two-dimensional invariant torus 
is born from the limit cycle. The third Andronov-Hopf bifurcation occurs when a three-dimensional invariant 
torus is born from the two-dimensional one (filled by quasiperiodic orbits) . The so-called Landau- Hopf scenario 
of the onset of turbulence envisioned a chain of further Andronov-Hopf bifurcations which would lead to a 
creation of an invariant torus of a sufficiently high dimension, i.e. to a quasiperiodic regime with a high number 
of rationally independent frequencies (see more discussion in |23j). However, as Ruelle and Takens pointed out 
in PP, the dynamics on the invariant torus is not necessarily quasiperiodic: at the moment the torus is born, the 
system can be perturbed in such a way that every orbit on the torus will become periodic, hence the Poincare 
map will be the identity, hence, as it follows from our results above, a further small perturbation may lead to 
a chaotic dynamics for a flow on the torus of dimension 3 or higher. 

2 Universal maps in the Newhouse domain 

In general, it follows from Theorem [5] that every time we have a periodic orbit for which the corresponding 
Poincare map is, locally, identity: 

or coincides with identity up to flat (i.e. sufficiently high order) terms: 

x = x + o(||a;|| r ), 

a C"*-small perturbation of the system can make the Poincare map universal, i.e. bifurcations of this orbit 
can produce dynamics as complicated as it only possible for the given dimension of the phase space. Thus, 
arbitrarily complicated dynamical phenomena can be uncovered by the study of bifurcations of periodic orbits 
alone. 

In order to show how powerful this observation can be, let us apply it to the analysis of the dynamics of 
two-dimensional diffeomorphisms from the Newhouse domain. In the space of C r -smooth dynamical systems 
on any smooth manifold, Newhouse domain is the interior of the closure of the set of systems that have a 
homoclinic tangency (a tangency between the stable and unstable manifolds of a saddle periodic orbit; for a 
saddle periodic orbit to exist, the dimension of the phase space should be at least 3 in the case of continuous 
time and at least 2 for discrete dynamical systems). A non-trivial fact [TOJ [TTJ [5H [55] is that the Newhouse 
domain at r > 2 is always non-empty and adjoins to every system with a homoclinic tangency. Importantly, 
most of known global bifurcations which lead to the emergence of chaotic dynamics or happen within the class 
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of systems with complex (chaotic) behavior are accompanied by a creation of homoclinic tangencies. Therefore, 
Newhouse regions in the space of parameters can be detected for virtually every dynamical model with chaos 
(see more discussion in H7J [31 [7] ) . In my opinion, the ubiquitous presence of homoclinic tangencies in the 
dynamical models of a natural origin makes the study of maps from the Newhouse domain the most important 
problem of chaotic dynamics. It should also be mentioned that a commonly shared believe (with no hope to 
prove yet) is that the space of two-dimensional C r -difieomorphisms with r > 2 is the closure of the union of 
just three open sets: Morse-Smale systems (i.e. systems with simple dynamics), axiom A systems, and the 
Newhouse domain, i.e. unless a two-dimensional map with a chaotic dynamics is uniformly hyperbolic, it most 
probably lies in the Newhouse domain. 

Typically, a map from the Newhouse domain possesses an invariant basic hyperbolic set which is wild - the 
term meaning that for the map itself, and for every C r -close map (r > 2), there exists a pair of orbits within the 
hyperbolic set such that the unstable manifold of one orbit has a quadratic tangency with the stable manifold 
of the other [TQl [IT]. For a fixed pair of orbits, the corresponding tangency is a codimension-1 bifurcation 
phenomenon, so it can always be removed by a small smooth perturbation; the wildness, nevertheless, means 
that once the original tangency is removed a new tangency always appears, corresponding to some other pair 
of orbits from the same hyperbolic set. 

A wild hyperbolic set of a C r -diffeomorphism of a two-dimensional smooth manifold will be called ultimately 
wild if it contains a pair of periodic orbits such that the saddle value at one periodic orbit is less than 1 and 
at the other periodic orbit it is greater than 1. The open subset of the Newhouse domain which corresponds 
to maps with ultimately wild hyperbolic sets will be called the absolute Newhouse domain. The saddle value 
is, by definition, the absolute value of the product of multipliers of the periodic orbit. If Q is a period q 
point of a map / (i.e. f q Q = Q), its multipliers are the eigenvalues of the matrix of derivatives of f q at the 
point Q. Therefore, if the saddle value is greater than 1 in the absolute value, then the map / expands area 
near Q, and it is area-contracting near Q if the saddle value is less than 1. So, no map from the absolute 
Newhouse domain is area-contracting, nor area-expanding. Moreover, the persistent tangencies between the 
stable and unstable manifolds of the wild hyperbolic set mean that none of these maps is uniformly hyperbolic, 
nor uniformly partially-hyperbolic. Thus, there is no obvious restrictions on the dynamics of two-dimensional 
diffeomorphisms from the absolute Newhouse domain, and the following theorem is, therefore, in agreement 
with the general "guiding principle" formulated in Section [TJ 

Theorem 3. For every r > 2, a C r -generic diffeomorphism from the absolute Newhouse domain is C r - 
universal. 

Proof. Fix any r > 2 and let a C-diffeomorphism of a smooth two-dimensional manifold possess an 
ultimately wild basic hyperbolic set A. Let P and Q be two saddle periodic points in A (f p P = P and 
f q Q — Q) with the saddle values a p and o Q such that o p < 1 and a Q > 1. As the periodic points P and Q are 
hyperbolic, they are preserved at all sufficiently small perturbations. Recall also that the unstable manifold of 
every point in the basic hyperbolic set has a transverse intersection with the stable manifold of every other orbit 
in this set. Therefore, the invariant unstable manifold W U (P) intersects the invariant stable manifold W S (Q) 
transversely at the points of some heteroclinic orbit Tpq. By the transversality, this orbit is preserved for all 
maps sufficiently close to /. According to [6], given every m > 1, in any neighborhood of / in the C r -topology, 
there is a C^-diffeomoprphism for which the unstable manifold W U (Q) has a tangency of order m to the stable 
manifold W 8 {P) at the points of some heteroclinic orbit Tqp (Fig. 3). This is a non-trivial statement: while the 
possibility to obtain a quadratic tangency (i.e. m = 1) by an arbitrarily small perturbation follows immediately 
from the wildness of A and from the fact that the stable and unstable manifold of any given periodic orbit in 
the basic hyperbolic set are dense within the union of stable and, respectively, unstable manifolds of all orbits 
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in this set [TU1 E], the tangencies of higher orders are due to the existence of moduli of local f2-conjugacy 

lump. 

Let / G C°° be a C r -close to / diffeomorphism for which the above described heteroclinic cycle exists. This 
cycle is a closed set that consists of 4 orbits: two periodic orbits (the orbits of P and Q) with a p < 1 and 
<j Q > 1, a transverse heteroclinic Tpq and the orbit Tqp of heteroclinic tangency of a sufficiently high order 
to. We will call an open region filled by periodic orbits of the same period a periodic spot. In Section [5] (see 
Lemma [3] and Remark after it) we show that 

for any neighborhood U of the heteroclinic cycle, arbitrarily C r -close to f there is a diffeomorphism f which 
has a periodic spot whose all iterations lie in U . 

For every orbit in the periodic spot the corresponding Poincare map (the map f k where k is the period of 
the spot) is, locally, identity x = x. Hence, as we mentioned in the beginning of this Section, it follows from 
Theorem [2] that in any C"*-neighborhood of / there exist C"*-universal maps. 

We have just shown (modulo Lemma [3] which is proved in Section [5]) that universal maps are dense in 
the absolute Newhouse domain. This immediately implies the genericity of the universal maps. Indeed, given 
an orientation-preserving C r -diffeomorphism g of the unit disk into R 2 , denote as V(g, 6) the set of all C r - 
diffcomorphisms from the absolute Newhouse domain whose dynamical conjugacy class intersects the open 
^-neighborhood of g (i.e. whose certain renormalized iteration is at a C r -distance smaller than S from g). This 
set, by definition, contains all universal maps — hence, it is dense in the absolute Newhouse domain. This set 
is also open by definition. Take a countable sequence of maps gi which is dense in the space of orientation- 
preserving C '-diffeomorphisms of the unit disk into R 2 , and a sequence Sj of positive reals converging to zero. 
By construction, the countable intersection V(gi, Sj) is a residual subset of the absolute Newhouse domain, 
and every map that belongs to this set is universal. □ 

In essence, this theorem gives an exhaustive characterization of the richness of dynamical behavior in the 
absolute Newhouse domain: every two-dimensional dynamics can be approximated by iterations of any generic 
map from this domain. In a simpler case of the Newhouse domain in the space of area-preserving maps, a 
similar statement is contained in [7]: iterations of a generic area-preserving map from the Newhouse domain 
approximate all symplectic dynamics in a two-dimensional disc. For area- contracting maps, it follows from 29 
that the closure of the dynamical conjugacy class of a generic map from the Newhouse domain contains all 
one- dimensional maps (we cannot have truly two-dimensional dynamics here, as the areas are contracted). As 
we mentioned, any two-dimensional map which is not in the class of maps with a uniformly-hyperbolic structure, 
nor on the boundary of this class, falls, hypothetically, in one of the three types of the Newhouse domain: the 
first is filled by area-contracting maps, the second by area-expanding maps (i.e. the maps inverse to the maps of 
the first type), and the third is the absolute Newhouse domain. We see that our Theorem [3] somehow completes 
the description of two-dimensional dynamics. 

The three types of Newhouse domains of two-dimensional maps were introduced in [30) . It has been known 
since [TU] that a generic area-contracting map from the Newhouse domain has an infinite set of stable periodic 
orbits, and the closure of this set contains a (wild) hyperbolic set. The latter fact is especially important: 
while chaotic dynamics is usually associated with hyperbolic sets, i.e. with saddle orbits, the Newhouse result 
shows that stable periodic motions can imitate chaos arbitrarily well, and they indeed do it generically. In 
[30] . for the Newhouse domain of the third type, it was shown that a generic map has both an infinite set of 
stable periodic orbits and an infinite set of repelling periodic orbits, moreover the intersection of these sets is 
non-empty and contains an ultimately-wild hyperbolic set. Thus, not only the Newhouse phenomenon holds, 
we also have a new effect here: a generic inseparability of attractors from repellers. In fact, these attractors and 
repellers can be more complicated than just periodic orbits: in |31j . the coexistence of infinitely many closed 
invariant curves was established for two-dimensional maps from the Newhouse domains of the third type. Our 
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Theorem [3] strengthens these observations: it implies the coexistence of infinitely many hyperbolic attractors 
and repellers for a generic map from the absolute Newhouse domain (see remarks to Theorem [2}. As we obtain 
the hyperbolic attractors and repellers by a perturbation of periodic spots, and the periodic spots are found 
in an arbitrarily small neighborhood of any heteroclinic cycle of the type we consider in Theorem [31 it follows 
that the closures of the set of hyperbolic attractors and of the set of hyperbolic repellers that we construct here 
contain any transverse heteroclinic orbit connecting the points P and Q. Such heteroclinic orbits are dense in 
the basic hyperbolic set A. Thus, it follows from the proof of the theorem that 

a C r -generic two-dimensional map from the absolute Newhouse domain has an infinite set of uniformly-hyperbolic 
attractors and an infinite set of uniformly-hyperbolic repellers, and the intersection of the closures of these sets 
contains a non-trivial hyperbolic set. 

In what follows we prove Theorem 1 and finish the proof of Theorem [3J 

3 An approximation theorem. 

Let F be an orientation-preserving C-diffeomorphism (r > 3) which maps the ball B n : {X)™=i x 1 — 1} 
R n , Without loss of generality we may assume that F is extended onto the whole R n , i.e. it becomes a C r - 
diffeomorphism R n — > R n , and it is identical (i.e. F(x) = x) at ||x|| sufficiently large; such extension is always 
possible. Let if be a constant such that 

(3.1) 

xeR n J(X) 

where J{x) is the Jacobian of F. Denote i?" := {x n > 0}. 

Lemma 1 . There exists a pair of volume-preserving, orientation-preserving C r ~ 2 -diffeomorphisms $i : i?" — > 
Rl and $ 2 : R n -> R n such that 

F = $2°*2 ^i°tii (3.2) 

where 

*j := (xi,...,x n -i,x n ) H> (xi,...,x n -i,ipj(x n )) (J = 1,2), (3.3) 

with 

V>iO„) = e Kx '\ ^ 2 {x n ) = Inav (3.4) 

Proof. We need to construct a volume-preserving diffcomorphism $i : (xi, . . . , x„ > 0) H> (xi, . . . ,x n > 0) 
in such a way that 

det— ^ 2 o$i o* x (x) = J(x) (3.5) 
ox 

(then the Jacobian of $2 = F o (^2 o $! o ^i) -1 will be equal to 1 automatically). By (|3.3p . (13.4l) . condition 
(13. 5p is equivalent to 

— if \ Kx n 

x n (pyzi i • ■ • i x n ) — j , 1 1 \ ' 

J yXi , . . . , X n — i , m x n J 

It follows from (|3.ip that 

d(f>/dx n > (3.6) 
everywhere. Moreover, as F is the identity map outside a bounded region of R n , we have that 

4>(x) = Kx n (3.7) 
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outside a compact subregion of i?" . Therefore, every trajectory of the vector held 

d(j) dcf) 

x 3 =0 (j<n-2), = - — , x n = ~ (3.8) 

ax n OX n -l 

extends for all x n -i G (—00, +00), and the time t(x) that the trajectory of the point x needs to get to x n -\ — 
is a C r_2 -function of x, well defined everywhere in i?" . Moreover, as it follows from (j3.7[) . (|3.8[) 

t(x) = -—x n -i + t q (x), (3.9) 

A 

where to is a uniformly bounded function, vanishing identically at x n close to zero and at sufficiently large 
x n . Thus, for every fixed values of Xj (j < n — 2), given any C £ (— 00, +00) the level line t(x) — C in the 
x n ) -plane coincides with the straight line x n -\ = —KC at x n close to zero and at sufficiently large x n . 
Every such level line is a connected set (as it is the image of the line {x n _i — 0, x n > 0} by the time-(— C) map 
of the flow of (I3.8p ). Thus, as x n runs from to +00, the value of <f> on this line runs all the values from to 
+00 (see (|3.7[) ). It follows that the map i?" — > i?" defined by 

x i = x 3 (j<n-2), x n -i = —T(x), x n = 4>(x) (3.10) 

is surjective. 
By O, 

9t d(f> _ dr_ d(f> = _ J 
dx n _i dx n dx n dx n -i 

It follows that for every fixed values of Xj (j < n — 2), the function ip changes monotonically along every level 
line of r, which implies the injectivity of map (|3.10|) . Thus, map (|3.10[) is a C r ~ 2 -diffeomorphism i?" — )• i?™. 
By (|3.1ip . it is volume-preserving and orientation-preserving, i.e. it is the sought map $1. □ 
The maps x 1-^ x of the following form: 

Xl = X 2 , ■ ■ ■ , Xn-l = %n, X n ^ (-l) n+1 X 1 + h(x2,--.,X„) (3-12) 

(note no dependence on xi in h), will be called Henon-like volume-preserving maps. Note that such maps are 
always one-to-one, and the inverse map is also Henon-like. 

Theorem 4. Every orientation-preserving C ' -diffeomorphism F : B n — > R n can be arbitrarily closely approxi- 
mated, in the C-norm on B n , by a map of the following form: 

H 2q2 o...off 2 iof 2 o H Ul o . . . o Hu o (3.13) 

where the maps are given by t3.3\) . and Hj s (j — l,2;s = l,...qj) are certain polynomial Henon-like 
volume-preserving maps. 

Proof. First, take a C r+2 -diffeomorphism F which approximates F sufficiently closely in C r . For the map 
F construct the decomposition F = $ 2 ° ^2 ° $1 ° *f i given by Lemma [TJ all the maps in the decomposition 
are at least C r . The map $i can be extended onto x n < by the rule x n — Kx n , x n -i — x-n-i/K (see 
(|3.I0p . (|3.7p . (|3.9p ). so it becomes a volume-preserving, orientation-preserving C-diffeomorphism R n — > R n . 
Then the theorem follows immediately from Lemma [2] below. □ 
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Lemma 2 . Every volume-preserving, orientation-preserving C r -diffeomorphism <I> : R n — >• R n can be arbitrarily 
closely approximated, in the C -norm on any given compact, by a composition of polynomial Henon-like volume- 
preserving maps. 

Proof. At n — 2 this result is immediately given by Theorem 1 in [3], so we proceed to the case n > 3. 
It is well known that $ can be imbedded in a smooth in t family Tt of volume-preserving C-diffeomorphisms 
R n — > R n such that = id and T\ = $. The derivative j^J-t defines a divergence-free vector field X(t,x), 
i.e. the diffeomorphism J- t is the time-i shift by the flow generated by the field X. One can approximate 
X arbitrarily closely on any given compact by a C°°-smooth divergence-free vector field which is defined and 
bounded for all (x,t) € R n x [0, 1]. Therefore, it is enough to prove the lemma only for those $ which can be 
obtained as the time-1 shift by the flow generated by such a vector field, i.e. we may assume that X € C£° with 
no loss of generality. 

Let T T t = J-t+T o , i.e. it is the shift by the flow of X from the time t to t + r. This map is 0(r)-close 
to identity, in the C r -norm on any compact subset of R n x [0, 1]. By construction, given any arbitrarily large 
integer N, 

$ = T t,(n-i)t 0...0 T T . mT o . . . o T t .q (3-14) 

where r = l/N, and m = 0, . . . , N — 1. 

Note that the vector field X admits the following representation: 

n-i 

X = J2 xii) ( 3 - 15 ) 

i=l 

where X^ is a C°°-smooth divergence-free vector field such that 

Xj=0 at (3.16) 

Indeed, if we write the field X as 

X\ (x, t), i 1,...,ti, 
where -r— ^- = (the zero divergence condition), then the fields are defined as 

1 OX% 

with 

»7i=fi> Vi = ^i-Ci-i (i = 2,...,n- 1), 

/•Xi+l g 

(i = - ——ri l (xi,...,Xi,s,x l+ 2,...,x n ,t)ds (i = 1, . . . ,n - 2), ( n _i = £„. 

By construction, the fields X^, . . . ,J>f("~ 2 ) are divergence-free, and X^ 1 ' = X - - ... - X^" 2 ), so 
X(n-i) j g a j g0 divergence-free, as X is. 
It follows from (|3.15p that 

Tn = Tir l) o . . . o T« + 0(r 2 ), (3.17) 

where is the shift by the flow generated by the vector field jW. Recall that the maps T T ^ T in (|3 . 14[) are 
0(l/JV)-close to identity. Therefore, it follows from (IXT71) . (|3~TH) that 

* = ^ 1)T ° • ■ • ° Tl%_ 1)T 0...0 o . . . o T±% T o . . . c o . . . o T$ + O(r), (3.18) 



10 



uniformly on compacta. 

As t can be taken arbitrarily small, it follows that in order to prove the lemma, it is enough to prove that 

every of the maps T^f in the right-hand side of (ETTg|) can be approximated arbitrarily well by a composition 

(i) 

of Henon-like volume-preserving maps. The maps T^ t are volume-preserving and satisfy 

Xj — xj at j ^ i, i + 1 (3.19) 

(see (|3.16[1 ). Therefore, if we denote 

Xi=p(x), x i+1 = q(x), (3.20) 

then 

® M =1. (3.21) 



d(xi,x i+1 ) 

Thus, we can view (|3.20[) as an (n — 2)-parameter family of symplectic two-dimensional maps (2^,2^+1) 1— > 
(xi,x i+ i) parametrized by (xi, . . . , x l -i,x i+2 , ■ ■ ■ ,x n ). 

According to 3], every finite-parameter family of symplectic maps can be approximated (on any compact) 
by a composition of families of Henon-like maps, i.e., in our case, maps of the form 

Xi — Xj-j-i , X^+i X. L -\- /l(x^_|_i,Xi, . . . , , Xi^-2 j ■ • ■ ; *^n) ■ 

It follows that every map of the form (|3.19p . (l3.20[) . (l3.21j) can be approximated arbitrarily closely by a compo- 
sition of the maps of the form 

Xj = Xj at j ^ 1, 

Xi=x i+1 , (3.22) 

Xi+l Xi -\- h(xi-\-i , Xi , . . . , Xi—i , Xi-^2 ) ■ • ■ ; X71 ) ■ 

It just remains to note that every map of form (|3.22|) is a composition of volume-preserving Henon-like maps; 
namely, it equals to 

S 11 - 1 - 1 o H o S o Q' 1 - 1 o S l+1 , 

where 

S := (xi, . . . ,x n ) h-> (x 2 , . . . ,x n , (-l) n+1 xi), 

n 

Q := (xi,...,x n ) H' (x 2 , ■ ■ ■ ,x n , ^2(-l) n+j Xj), 

3=1 

H . X2, • • • j %n — 1 — $711 3-"n — 

n-1 

= " Xn + h(x n] x n _ t+1 , . . (-1)" +1 X 2 , . . . , (-1)™ +1 .T„_ 4 )}- 

3=1 

End of the proof. □ 
Remark. Consider the map 

$o := (a;i,...,a;„_2,a; n _i,Xn) >-> (x%, . . . ,x n ^ 2 ,x n , -x n -%). (3.23) 
This is an orientation- and volume- preserving diffeomorphism of i?™, therefore we may rewrite (I3.2[) as follows: 

F = $ool 2 oit 2 o$ (l o$ 1 o vj/^ 
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where $12 are orientation-preserving, volume-preserving C'~ 2 -diffeomorphisms ($j = $q 1 o we assume 
that $1 is extended onto the whole of R n , like in Theorem|4]). Now, by Lemma[2l we obtain the following, more 
convenient for us, analog of Theorem @] 

the map F can be arbitrarily closely approximated by a map of the following form: 

$o ° H 2q2 o . . . o H 21 o 1< 2 o $ o H lqi o . . . o H u a (3.24) 
with polynomial Henon-like volume-preserving maps Hj S . 

4 Proof of Theorem 1. 

Given a diffcomorphism F : B n — > R n we will take its sufficiently close approximation F in the form of (|3.24p . 
Then we will construct a close to identity map (which we denote Yg below) whose some renormalized iteration 
is a close (as close as we want) approximation to F. 

The map Yg is a small perturbation (as small as we want, in the C r -norm for any aforehand given r) of the 
time-<5 map Yg of a certain C°° flow Y in R n ; the constant S can be chosen as small as we need (i.e. Yg is indeed 
a small perturbation of the identity). In our construction, the vector field of Y vanishes identically outside some 
ball D that does not depend on the choice of the approximation F, it does not depend on S either. The small 
perturbations which we will apply to Yg will also be localized in D. Hence, our close to identity maps Yg are 
all equal to identity outside the same ball D. Since their renormilzed iterations approximate F arbitrarily well, 
this gives us Theorem 1, as we explained it in the remark after the theorem. 

We define the flow Y by means of the following procedure: we give explicit formulas for the vector field inside 
certain blocks U\±, U 2 ±, Via described below, while between the blocks we specify only the transition time from 
the boundary of one block to another and the corresponding Poincare map. The existence of a C°° flow with 
arbitrary (of class C°°) transition times and orientation-preserving Poincare maps between block boundaries is 
a routine fact (at least for the given geometry of the blocks, see Fig. 2). 

The idea of the construction is as follows. Given a diffcomorphism F : B n — > R n , we may always approximate 
it by a C°°-diffeomorphism arbitrarily well, so we assume that F g C°° from the beginning. Let $12 and $12 
be the maps defined by decomposition (|3.2p of F. We define the vector field inside the blocks Ui±, U 2 ± in such 
a way that a kind of saddle- node bifurcation is created inside each of the blocks (see (|4.1.1l) . (|4.2.ip ; we build 
very degenerate saddle- nodes in order to make formulas for the time-i map simpler - see (|4.1.3[) . (|4.2.4[1 ). The 
Poincare maps from the boundary E°™' of Uj+ to the boundary E™_ of Uj- (j = 1,2) are chosen in such a 
way that the resulting flow map from entering Uj+ till exiting Uj- is linearly conjugate to the map ^ j (see 
(I4T51) . (I4X5|) . (1331) .(1X41)1. 

We make the flow in Vx$ volume-preserving and linear (see (|4.1.6[) . (|4.2.2j) . (|4.2.3jl ). Moreover, we put saddle 
equilibria into Via- We define the Poincare map between the boundaries 11°"* and ITj™ of Vj (see Fig. 2) in such 
a way that a homoclinic loop to the saddle is created. For the time-(5 map Yg of the flow Y the saddle equilibrium 
is a saddle fixed point, and the homoclinic loop is a continuous family of homoclinic orbits. We perturb the 
map Yg in such a way that this family splits into a finite set of orbits of homoclinic tangency of sufficiently 
high orders and unfold these tangencies then. The exact form of the perturbation (see (|4.1.17p . (l4.2.16p ) may be 
chosen such that the iteration of the perturbed map Yg that corresponds to one round near the homoclinic loop 
is a close approximation (in appropriately scaled coordinates) to any given polynomial conservative Henon-like 
map (see (|4.1.22p . (|4.2.22p ). Hence, a multi-round iteration of Yg can be made arbitrarily close to a superposition 
of a finite number of any given polynomial Henon-like maps. By Lemma [2] such superpositions approximate any 
given volume-preserving maps, e.g. maps $j from the decomposition (|3.2p . In this way an iteration of the map 
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Yg that, after a large number of rounds near homoclinic orbits, takes points entering Vj to the points entering 
Ua-j is made as close as we want to the map §j (j — 1,2), in some rescaled coordinates. 

The parameters of the flow and perturbations are chosen in such a way that the coordinate scalings we make 
in the blocks Ut,2 and Vt,2 match each other. Thus, by construction, the (renormalized) iteration of Yg that 
corresponds to passing from the entrance to Ut+ through U\- into Vt, then to many rounds near the homoclinic 
loop, then to exiting V± and entering U2+, passing through U2- to V2, a number of near- homoclinic rounds and 
return to Ui+, is a close approximation to $2 ^2 $1 ^1, i-e. to the original diffeomorphism F indeed. 

In the two-dimensional case, the key fact that every area-preserving diffeomorphism can be approximated 
by some multi- round iteration of a map with a homoclinic tangency was proven in [6] . To deal with dimensions 
higher than 2, we construct a very degenerate homoclinic tangency (in terms of [20) . both critical and Lyapunov 
dimensions for this tangency are equal to n - this is necessary in order not to lose dimension at the rescaling 
of the first- return map). We do not undertake an analysis of the corresponding bifurcation; instead, we make 
explicit computations of the rescaled return maps for our particular example only. 



4.1 Two-dimensional case 

To make computations more transparent we start with the case n = 2. Let $1,2 and ^1,2 be the maps defined 
by (|3.2j) . Let I\± and fo± be intervals of values of X2 such that S I\+ at (xt, X2) € B 2 , X2 € It- at (xt, X2) € 
^i{B 2 ), x 2 £ h+ at (xt,x 2 ) e $io*i(B 2 ) and x 2 € h- at (xt,x 2 ) € * 2 ° $1 ° ^>i{B 2 ). Let R be such that all 
the intervals Ij± lie within {\x2\ < R}- Choose numbers at+ = ai_+3 = 6i + 6 = a2+ + 9 = 02- + 12 = &2 + 15. 
Let the vector field of Y in the regions Uj a : {\xt — a Jcr || < 1, \x%\ < R}, j = 1, 2, a = ±1, be equal to 

Xt = - (1 - MiX 1 ~ £Ol ~ O-ja)), u l x x 

±2 = <7Jj X 2 £(xt ~ djo-), 

where ^1,2 > are small (see (|4.1.3ip ~). 71,2 € [0, 1] (see (|4. 1.321) ). and £ is a C°° function such that 

0<£<1, £(0) = 1, t(z) = at \z\>± (4.1.2) 

As ii < in Uj a , every orbit of Y that starts in Uj a near x\ = a,j a + 1 must come in the vicinity of 
x\ = a.j a — 1 as time grows. For the corresponding time-i map, we have 

xt(t)=xt(0)-t+^p{^), x 2 (t) = e^ a ^x 2 (0), (4.1.3) 



where (see (|4.1.1|l ) 



l W J x „_ 1 (t)-a^ J-l/2 jU-Kl-ZiXl-^z))' 

(4.1.4) 



- 2 (j^/ 2 /i+C i_ A o( 1 _ (£(js)) 
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Note that a(/x) is positive, independent of 2i(0) and t (because we assume that the integration limits (xi(0)— dja) 
and (xt(t) — a,j a ) are close to 1 in the absolute value, i.e. they fall in the region where £(z) = 0; see l|4. 1.2ft ). 
and both a{fj) and /3(/it) tend to infinity as n — > +0 (the integrals diverge at /i = because £(0) = 1). 
Denote £}" := { Xl = a J+ + 1, |ar 3 | < 1}, E£* := {<Bi - a i+ - 1, |z 3 | < i?}, Ejl := {a* = Oj- + 1, |a? a | < 
:= {xi = <ij_ — 1, 1 353 1 < 1}- Every orbit of Y that intersects ££j_ at £2 sufficiently small leaves f7j + by 
crossing £°" 4 , and the orbits that intersect leave £/j_ by crossing £°™' (see (|4. 1.3ft ). Define the vector field 
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Figure 2: An illustration to the proof of Theorem 1. 
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Y in the region between E°"' and Ej™ in such a way that the orbits starting in E°_"' reach E 1 ™ at time 1, and 
the corresponding Poincare map E°"' — > EJI is 

.T 2 !->■ 1pj(x 2 ), 

where ipi(x 2 ) = e^ 12 at x 2 S ii+ and ip 2 (x 2 ) = hix 2 at x 2 G 7 2 + (see (I3.4[l ). Then, the flow takes the points 
from the vicinity of x± = Oj+ + 1 in Uj+ into the vicinity of x\ = aj- — 1 in Uj—. By (14.1.31) . the corresponding 
time-i map Sjt is 

*i(t) = x 1 (0)-t + P(p j ), x 2 (t) = e -Ti«(w) V,- (e^ Q ^x 2 (0)) . (4.1.5) 

In the regions Vj : {|xi — bj\ < l,\x 2 \ < 1} (j = 1, 2) we put the vector field of Y to be equal to 

±1 = — (xi — bj), x 2 =x 2 . (4.1.6) 

Thus, in Vj, the point Oj : {x\ — bj,x 2 — 0} is a linear saddle. Its local stable manifold Wf is x 2 = 0, and the 
local unstable manifold VP" is x\ = bj. The time-i map Lj t within Vj is given by 

xi(t) = b j + e- t (xi(0)-bj), x 2 (t) = e'x 2 (0). (4.1.7) 

Let us define Y in the region between E°"' and IE™ := {x\ = bj + 1, \x 2 \ < 1} in such a way that all 
the orbits starting in a small neighborhood of x 2 = in E°"' intersect ITj™ at time 1, and the corresponding 
Poincare map is the identity: 

x 2 i-> x 2 - 

Then the time-i map Qj t from a small neighborhood of xi = Oj_ — 1, x 2 = in £/j_ into a small neighborhood 
of ^! = bj + 1 in is given by 

ari(t) - bj = e -(*-*i(o)-a+ai-) ) ^(j) = e *- a ' l ( o >- a+, y-a:a(0) (4.1.8) 

In order to see this, we recall that the vector field in Uj- near x\ = aj— — 1 is given by 

xi = —1, ±2 = 

(see (|4.1.1|l ). Therefore, the term xi(0) + 2 — <ij_ in (|4. 1 .8[) is the time the orbit spends in order to get from 

x(o) to ny; ). 

Every orbit that enters Vj at x 2 > leaves Vj by crossing the cross-section 11°"' := {x 2 = 1, — fej| < 1}, 
and every orbit that enters Vj at x 2 < leaves it by crossing the cross-section 11°"' := {x 2 = —1, \x% —bj\ < 1}. 
We assume that the orbits that start at Hj+ close to the point VP" n 11°"* = (bj, 1) return to Wj at time 1 and 
cross ITj™ := {x\ =bj — 1, \x 2 \ < 1}; we also assume that the corresponding Poincare map x\ ^ x 2 is given by 

x 2 = —(xi- bj) 

(the minus sign stands to ensure the orientability of the flow map). It follows that the time-t map Tj t from a 
small neighborhood of (~l 11°"' in Vj into a small neighborhood of X\ = bj — 1 in Vj is given by 

Xl (t) = bj - e-^-V /a 2 (0), ar 2 (*) = -e t - 1 a:2(0) a (xi(Q)-6 i ). (4.1.9) 
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For the orbits that start at 11°"* close to the point W" n 11°"' = (bj, -1) we assume that they cross £|V. j+ 
at time 1, and the corresponding Poincare map is given by x 2 — —{x\ — bj). Thus (see (|4.1.ip . (14.1.2p , (|4.1.6[) ). 
the time-i map Gjt from a small neighborhood of Wf 011°"' in Vj into a small neighborhood of X\ — a?>-j,+ + 1 
in U 3 -j,+ is 

xi(t) = a 3 - jt+ + 2 - t - In 1x2(0)1, z 2 (t) = -(zi(Q) - b 3 )\x 2 (0)\. (4.1.10) 

Every C°° flow F, which satishes (|4.1.5p . (|4Tl~S)) . (|4.1.7|) . (j4TO1) . (j4.1.10l) . is good for our purposes. We may 
therefore assume that the vector field of Y is identically zero outside some sufficiently large ball D. For small 5, 
the time-5 map Y$ of the flow is 0(5)-close to identity in the C"*-norm, for any given r. It also equals to identity 
outside D. Let us fix a certain r, and take a sufficiently small S (for convenience, we assume that N := S^ 1 
is an integer). Below we construct an arbitrarily small (in the C r -norm), localized in D perturbation of Y$ as 
follows. 

For the given diffcomorphism P, take its sufficiently close approximation in the form of (I3.24[) ; in the two- 
dimensional case the map $o is given by 

$o := (xi,x 2 ) H> (x 2 ,-xi). (4.1.11) 

As there is only a finite number (q\ + q 2 ) of the polynomial maps Hj S in (|3 . 24[> . one can find some finite d > 1 
(common for all Hj S , s — 1,. . . ,qj,j — 1,2) such that the maps Hj S are written as follows: 



xi = #2, x 2 = —x\ + ^ h jsv x 2 . (4.1.12) 



0<i/<d 

In the segment I° ut :— {e~ s < x 2 < 1} of Wf (j = 1, 2), we choose qj — 1 different points Mji, . . . , Mj iq ._i, 
and one point Mj qj £ W" will be chosen in the segment — e~ s > x 2 > — 1. Let Uj s denote the coordinate x 2 of 
M js (s = 1, . . . ,^). As iV<5 = 1 =flight time from n°f to II™ , near the segment I° ut the (N + l)-th iteration 
of the time-i5 map Ys is the map Tj,i+,5 from (|4. 1 . 9[) . i.e. it is given by 

x\ = bj — e~ s /x 2 , x 2 = —ex\(x\ — bj). (4.1.13) 

This map takes the segment I° ut onto the segment {bj—1 < X\ < bj—e~ s , x 2 = 0} £ W?. Let Pj, s +i = Tj^ + $Mj s 
(s = 1, . . . ,qj — 1), and let Pj\ be a point from {bj+e~° < x\ < bj + 1, x 2 = 0} £ WJ. We denote the coordinate 
xi of Pj.s+i as z js . By (I4.1.13[) . 

Zj, s +i = bj - e- s /u js . (4.1.14) 

We also assume 

Zjx = bj + e- s/2 . (4.1.15) 

We stress that Uj s and (6j — Zj s ) are bounded away from zero. 

Take a sufficiently large integer m and choose some points Pj s = (zj S ,z^ s ) and Mj s = (Uj- a ,Uj S ), sufficiently 
close to P, s and M,- s respectively (j = 1, 2; s = 1, . . . , cfy). We define the coordinates of Pj s and Mj s by the 
following rule: 

z ' js = e- m u js ), u' js = bj + e- m {z JS - bj) at s < qj - 1 

/ =e -{m+l j 8) I _ b . +e -(m+l 3 S) { _ b . ) (4.1.1b) 

where lj is an integer to be defined later (see (|4.1.24[) : note that IjS remains uniformly bounded for arbitrarily 
small 5). As m is assumed to be large and IjS is bounded, such defined points Pj s and Mj s are closed to 
respective points Pj S and Mj S indeed. It follows from (|4. 1 . 16|) that Mj s = Lj m Pj s at s < qj — 1, where Lj t is 
the map (|4. 1 .7j) . At s = q 3 we have M' jqj = Lj jm+ljS P'j q . . 
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Let us add to the map Y$ a small perturbation, which is localized in a small neighborhood of the points 
Yf 1 (Pj_ s+ i) (so outside these small neighborhoods Y$ remains unchanged). We require that these localized 
perturbations are such that in a sufficiently small neighborhood of Mj s — Y^ N (Yf 1 Pj^ s+ i) the map Y^ +1 
(where Y$ denotes the perturbed map) is given by the following perturbation of (14.1.13[) : 

xi = bj - e~ s /x 2 , 

x 2 = z' ha+l - e s xj(x 1 - Uj s ) + ^ s jsu(x2 - u js y , (4.1.17) 

0<v<d 

where Sj SL , are small coefficients to be determined later (see (|4.1.21|) V By (|4. 1 .16)) . z', ■ s+1 and (u'j s — bj) are 
small as well, so (|4.1.17[) is a small perturbation of (|4.1.13[) indeed. It is easy to see from (|4.1.17[) . (|4. 1.141) that 

At e = (hence at all small e) the map ° Lj m = Y^ nN+N+1 (where T stands for the perturbed map 

T) takes a small neighborhood of Pj s = (zj s ,z'j s ) into a small neighborhood of Pj a+ i = i z j,s+i, z'j s+1 ). By 
(|4.1.17p . (|4.1.7l) . (|4.1.16p . this map is written as 

x 1 = b j -e-( m+s >/x 2 , 

x 2 = z' jtS+l -e^ + ^x 2 2 (x 1 -z js )+ J2 e js ve um (x 2 -z' js y, (4.1.18) 

Q<v<d 

We choose some r](m) that tends to zero asm-> +oo and introduce rescaled coordinates (vi,v 2 ) near Pj s 
by the rule 

x\ = z js + C JS 77 vi, x 2 = z' js + — i] e' m v 2 , (4.1.19) 

where the independent of m positive coefficients Cj S are determined later (see (|4.1.20[) ; note that Cj S are 
bounded away from zero and infinity). Since 77 tends to zero as m — > +00, any bounded region of values of v 
corresponds to a small neighborhood of Pj s . 

After the rescaling, map (|4.1.18j) takes the following form (see (|4.1.16p . (|4.1.14p ): 

Cj, s +lVl = -JjT7\ v 2 + 0(f)), 

1 



' " U 'js C 3s 



-v 2 



s+l 



= -e s u%C js v 1 +0(r ] )+ £ Ejav e m rf- 1 C^, 



0<u<d 



As we see, by putting 



and 

the map Xj,i+<5 o Lj m near Pj s takes the form 



= (4 - L2o) 

s.« = '"'/ ; "c,.:. /I.- (4.1.21) 



vi = v 2 + 0( V ), 1 
v 2 = -vi + 0(77) + £o <iy <d h jsv v%, 



i.e. it can be made as close as we want to the map Hj S , provided m is taken large enough (recall that 77 — > as 
?7i — > +00). We take 77 tending to zero sufficiently slowly, so all £j su — > (see (|4. 1 .21[) ; recall that v < d where 
d is independent of m and 5). Thus, our perturbation to Yg is arbitrarily small indeed. 
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It follows that in the rescaled coordinates the map fj^i+a ° Ljrr^j = y*( mAr+Ar+1 ) from a small neighbor- 
hood of Pj! into a small neighborhood of P^ q . can be made as close as we want to the map 



H M . o . . . o Hn, 



provided m is large enough. The rescaled coordinates near P^ and Pj g . are given by formulas ()4.1.19[) . where 
the coefficients Cji > may be taken arbitrary, and the coefficients Cj qj are then recovered from the recursive 
formula (|4.1.20p . Further it is convenient to put 

Cji = e- s/2 . (4.1.23) 

Note that m does not enter (|4.1.20|) , (|4.1.23|) , hence Cj qj stay bounded away from zero and infinity asm-> +00. 

Let us also fix the choice of the integer lj that enters the definition of the points Pj q . and M'j q . (see (|4.1.16p ). 
Namely, we require that 

ln|u^|+lnC^ =ljS. (4.1.24) 
Recall that Uj q j is the coordinate x 2 of the point Mj qj , and it can be arbitrarily taken within the interval 

-e~ s > u jqj > -1. (4.1.25) 

As Cj q . does not depend on the choice oiuj q j (see (|4.1.20|) ). condition (|4.1.24|) uniquely defines both the integer 
lj and the value of Uj q . . The claimed before uniform boundedness of IjS follows from the uniform boundedness 

Now, from (|4.1.7|l . (j4.1.19|) . ((4.1.16[) . ()4. 1 .24|) we obtain that the map Lj tm +>.s from a small neighborhood of 
Pj qj into a small neighborhood of Mj q . is identity (i.e. v± = v±, v% — v%), provided the rescaled coordinates 
(fi,^) near are introduced as follows: 

x x = u' jq . + \u jqj T 1 77 e' m v u x 
X2 = Uj qj -u jqj T) v 2. 

Therefore, for the given choice of the coordinates, by taking to — > +00, the map Lj^+^s o ^Tj 1+ g o Lj m \ = 
f^ +mAr+ ^ (mAr+Ar+1) from a small neighborhood of Pj x into a small neighborhood of Mj 9 . can be made as close 
as we want to Hj qj o . . . o Hji (we have already proved the same for the map ^Tj t x + s o Lj m ^j , and the map 
Lj^m+ijS is identity in the chosen coordinates). 

Recall that, by construction, the point Y s N+1 M- q . lies in U^-j,^ in the region a 3 _ Ji+ + l — 5 < X\ < a 3 _ J:+ + l. 
We add to the map Yg a perturbation, localized near the point Y s N Mj q _, such that the corresponding map 
G-^i+i = Y S N+1 will have the following form near M' : 



xi = a,3-j, + + 1 - S - In \x 2 \, x 2 = -\x 2 \(xi — u'j q .); (4.1.27) 

since u'j q . — > bj as to — > +00 (see (|4.1.16[) ). these fomulas define a small perturbation of the map Gj^\+$ given 
by (|4.1T0|) indeed. 

Denote P$-j<a = Gj.i+sM'^.. By (|4.1.27p . the coordinates of ^3-^,0 are gi yen by {x\ — as-j :+ + 1 — S — 
In \uj qj I, x-2 = 0}. Introduce rescaled coordinates near the points P$-j (j = 1, 2) by the rule 

xi = a3-j,+ + 1 - S - In \uj qj I + 77 Vi, X2 = rje~ m v 2 . (4.1.28) 
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When coordinates are rescaled by rule (|4.1.26|) near Mj g , and by rule (|4.1.28[) near P^_j , map (|4. 1.271) takes 
the form (v\,v 2 ) *—¥ (v 2l — fi) + 0(rj), i.e. it becomes arbitrarily close to the map $o (see (|4.1.11j) ) as + oo. 
Thus, in the coordinates rescaled by rule (|4.1.19[) near P^ and by rule (|4.1.28l) near P^_j , the map Gj,i+s ° 

Lj^m+ijS ° {Tj^i+s o ij m ^ = y^ + ( 1+<? j)( mAr+A,+1 ) f r om a small neighborhood of P^ into a small neighborhood 

of P^-j o i can be made as close as we want to the map $o ° Hjqj ° ■ ■ - -ffji as m grows. 

Analogously, we take the point Mj : {x\ — aj- — 1 + (5/2, a; 2 = 0} € and perturb the map Y5 near 
Y s N Pj Q in such a way that the map Qj t i+g = near Mj will be given by 

Si - 6 rf = e-C^i-i+ai-), X2=e s - x ^- 1+a ^-x 2 + e- m u jl . (4.1.29) 

It is a small perturbation of the map Qj,i+s given by (|4. 1 .8[) . and it takes Mj to P^ (see (j4.1.15|l . (|4.1.16|l ). 
When we introduce rescaled variables near Mj by the rule 

x x =a,j- -1 + S/2 + r]Vi, x 2 =r 1 e~ m v 2 , (4.1.30) 

and near P'-^ by rule (|4.1.19|l . (|4.1.23|l . map (|4.1.29p will take the form v = v + 0(r/), i.e. it is close to the identity 

map. Thus, in the rescaled coordinates given by (I4.1.30[) , (|4.1.28[) , the map Gj^i+s Lj.m+^s ° ^T^i+s o L jm ^ o 

Qj,i+<5 = Yg 3+< " 1+q3 ^ mN+N+1 ' >+N+1 from a small neighborhood of Mj into a small neighborhood of Pjj.-q, is 
as close as we want to the map $0 Hjqj o . . . o if^ at m large enough, i.e. it is a close approximation of the 
map $ 3 . 

Let us now determine the form of the map Sjt '■ v 1— > v from a small neighborhood of Pj into a small 
neighborhood of Mj in the coordinates rescaled by the same rules (|4.1.28[) and (|4.1.30[) (note that when rescaling 
near Pj one should change j to 3 — j in the right-hand side of (I4.1.28[) ). We will choose the coefhcients fij in 
pXTj) such that 

fife) = fc - 5 + 3/25 + In |u 3 -j,g 3 _ j I (4.1.31) 

where k is some sufficiently large integer. Since j3 — > +00 as fj, — > +0 (see (|4.1.4p ). and In \u 3 -j m _ j | is bounded 
(see (|4.1.25[) ). equation (14.1. 3ip has a solution ^j(k) for every sufficiently large k, and ^j(fc) —> +0 as fc — ► +00. 
It follows that a(fij(k)) — > +00. Thus, for any sufficiently large m we can find G (0, 1] and large k such that 

T){m) e- m = e-T^CwW). (4.1.32) 

From (14.1.51) we immediately obtain that with this choice of fij, jj and k the map Sj^s = Y$ from a small 
neighborhood the point Pj into a small neighborhood of Mj takes the following form in the coordinates 
rescaled by rules, respectively, (|4. 1.281) (with (3 — j) changed to j) and (|4. 1 .30|) : 

vi = vi, v 2 = ipj(v 2 ). 

As wee see, the map Sj^s in the rescaled coordinates coincides with the map \t j for v from some open neigh- 
borhood of D (if j = 1) or of $1 o *x(D) (if j = 2). 
Summarizing, we obtain that the map 

G 2 ,i+s L 2: , n+ i 2 s o [t 2 s + s o L 2 m) o Q 2 ,i+5 S 2 ,ks° 

Gi^+j o L 1)TO+ ; l5 o ^Ti^i+i o L im ^ o Ql,l+<5 Sl,k« = 

_ y2k+l 1 +l 2 + (2+q 1 +q 2 )(mN+N+l)+ 2 (N+l) 
— 1 S 
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is a close approximation to the map $ 2 ° ^2 ° ° ^1 (i-c. to the original map F), provided Hj q . o . . . o Hj 1 
are sufficiently close approximations to <$>j (J = 1,2) and m is large enough. This completes the proof of the 
theorem in the two-dimensional case. 

4.2 Higher-dimensional case 

In the case of arbitrary n > 2, the construction follows the same line as in the two-dimensional case. As above, 
$1,2 and ^1,2 are the maps defined by (|3.2p . and I\± and I 2 ± are intervals of values of x n such that x n £ I\ + at 
x £ B n , x n e Ji_ at x £ ^i(B n ), x n £ I 2+ at x £ $1 o ^ 1 (B n ) and x n £ I 2 - at x £ * 2 o o ^(S"). A value 
i? is chosen such that all the intervals Ij± lie within {\x n \ < R}. Choose numbers a\ + = ai_ + 3 = bi + 6 = 
02+ + 9 = a 2 - + 12 = b 2 + 15. Define the regions Uj a ■ {\x n -i — a JCT || < 1, \x n \ < R, \xi\ < 1 (i < n — 2)}, and 
Vj : {|x n _i - bj\ < 1, < 1 (i^n- 1)}, j = 1, 2, ct = ±1. Let the vector field of a C°° flow Y in L/^ be 
equal to 

Xn-l = -fJ-j - (1 - - fC^n-l ~ a jcr))> 

(4.2.1) 

where m a > are small (see (|4.2.28|) . (|4.2.291> ). 7i± £ [0,1] (see (|4.2.29p . (|4.2.30p ). and the C°° function £ 
satisfies (|4. 1 .2j) . In the regions Vj we make Y equal to 

±i = -XiXi (i = 1, . . . ,7i - 2), i„_i = -A„_i(a;„_i - bj), x n = x n ; (4.2.2) 

here Xi > are such that 

A 2 = ... = A„_i = A, Ai=l-(n-2)A, (4.2.3) 

where the positive number A is specified below (see Q4.2.1ip ), 

As x n -\ < in Uj a , every orbit of Y that starts in Uj a near x n -\ = aj„ + 1 must come in the vicinity of 
x n —i — djc, — 1 as time grows. For the corresponding time-t map, we have 

x, l (t)=e^° a ^x t {0) (*<n-2), x n (t) = x n (0) - t + ^/3(p 3 ), (4.2.4) 

where the tending to infinity, as (1 — > +0, functions a(/j) and f3(/i) are defined by (|4.1.4|l . 

Denote E™ := {x„_i = a j+ + l,\x n \ < 1}, E°f := {x n ^ = a j+ - l,\x n \ < R}, E™ := {ir„_i = 
dj- + 1, |x n | < R}, E°"* := {x„_i = a,j- — 1, |x„| < 1} (we also assume that \xi\ < 1 for i < n — 2 on j]^ ou *)_ 
Every orbit of Y that intersects E2. at x n ,x%, . . . , x„_2 sufficiently small leaves Uj+ by crossing E™', and the 
orbits that intersect E*™_ leave J7,_ by crossing E°"' (see (|4.2.4p ). We define Y in the region between E°"' and 
E™. in such a way that the orbits starting in E™' reach E™ at time 1, and the corresponding Poincare map 
E°ut _^ jjot j g _ ^ x n _2, x«) i-> (xi, . . . , x n -2, ipj(x n )), where the functions tpj are defined by (|3.4p . Then, 
the flow takes the points from the vicinity of x n _i = aj+ + 1 in Uj+ into the vicinity of x n -\ = dj— — 1 in 
By (|4.2.4I) . the corresponding time-i map Sjt is 

x i (t) = e^+-^ a ^x i (0) (<<n-2), f . 

a: n _i(*) =a;»-i(0)-t + )8(Mi), a; B (t) = e-*i , »- a ^Vi(e 7 »+ a(M,) a: n (0)). 1 ' ' JJ 

In the region between E°"* and H"^_ := {x„_i = bj + 1, < 1}, we define Y in such a way that all the 
orbits starting in a small neighborhood of = in E°"* intersect IT*™ at time 1, and the 
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corresponding Poincare map is the identity: (xi, . . . , x„_ 2 , x„) H> (x%, . . . , x„_2, x„). Then the time-i map Qjt 
from a small neighborhood of x n _i = <ij_ — l,x„ = x\ = . . . = x„_2 = in Uj- into a small neighborhood of 
x n -% = bj + 1 in Vj is given by 



Xi(t) = e^-t'^-iM-^-fLfO) (i < n - 2), 

s„_i(t) - bj = e -M*-x»-i(o)-a+o J _) j ^ = e t-x„_ 1 (o)-2+a 3 _ Xn(0) 



(4.2.6) 



(see (gXU)). 

In Vjf (j = 1,2), the local stable manifold W? of the linear saddle equilibrium state Oj : {x„_i = bj,Xi = 
(i 7^ n — 1)} is x„ = 0, and the local unstable manifold Wj is x n -i = bj,x± = . . . = x n -2 = 0. The time-f 
map Lj t within Vj is given by 



Xi (t) = e- Xit Xi(Q) (i<n-2), , , 

x B _i(t) - 6, = e- A *(a;„_i(0) - 6,), x n (t) = e*x„(0). 1 ' ' ' 

Every orbit that enters Vj at x n > leaves Vj by crossing the cross-section 11°"* := {x n = 1, \x n -i — bj\ < 
1, |ar*| < 1 (« < n — 2)}, and every orbit that enters Vj at x n < leaves it by crossing the cross-section 
jjout ._ _ _^ |a> n _ 1 — 6j| < 1, |ari| < 1 (i < n— 2)}. We assume that the orbits that start at II "* close to 
the point W" nn°" 4 = {x\ = . . . = x„_ 2 = 0, x„_i = fej) return to Wj at time 1 and cross II*™ := {x n -i = bj — 
1, < l(i y^z n— 1)}; we also assume that the corresponding Poincare map (xi, . . . , x„_i) i-> (xi, . . . , x„_2, x„) 
is given by 

Xi = x i+ i («<n — 3), x n _ 2 = x„_i — 63, i„ = (-l)" +1 xi 

(the factor (— l) n+1 stands to ensure the orientability). It follows that the time-t map Tj t from a small neigh- 
borhood of Wj H I±°"* in Vj into a small neighborhood of x„_i = &j- — 1 in Vj is given by 

Xi (t) = e- A '( t - 1 )x„(0) A ^ 1 - A 'x l+1 (0) (i < n - 3), 

x„_ 2 (i) - e- A "- 2(t - 1) x n (0) A '- 1 - A - 2 (x n _i(0) - fy), , . 

x„_!(<) = 6, - e- A «- 1 (*- 1 )x„(0)- A »-S l4 -'' 8j 
x n (t) = (-l)"+ 1 e t - 1 x n (0) 1+Al x 1 (0). 

For the orbits that leave Vj by crossing 11°"*, we assume that the orbits that start at 11°"* close to the 
point n 11°"* = (xi = . . . = x„_2 = 0, x n _i = bj) cross £|™ - + at time 1, and the corresponding Poincare 
map (xi, . . . ,x„_i) 1 — ^ (xi, . . . ,x„_2,x„) is given by Xi = x, at i = 1, . . . , n — 2 and i„ = — (x„_i — Thus 
(see (|4.2.ip . (|4.1.2p . (|4.2.2|) h the time-f map Gjt from a small neig hborhood of Wf n II<™* in Vj into a small 
neighborhood of x„_i = 03_j 4. + 1 in is 

Xi(t) = |x n (0)| Ai Xi(0) (t<n-2), 

x„_i(f) = a 3 -j, + +2 - t -In |x n (0)|, (4.2.9) 
x„(t) = -|s n (0)| A »-*(a; B _i(0) -&J-). 

These conditions define the flow Y (we also assume that the vector field of Y is identically zero outside some 
sufficiently large ball D). Let us take a sufficiently small 5 such that N := S^ 1 is an integer, and proceed to 
the construction of a small (arbitrarily small in the C r -norm, with any chosen in advance r) perturbation of Yg, 
localized in D. 

Take a sufficiently close approximation of the given diffcomorphism F by the product (|3.24[) . There exists 
some finite d > 1, common for all Hj S (s = 1, . . . , qj,j = 1, 2) such that the polynomial maps Hj S in (13.24j) are 
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written as follows: 



Xi=x i+1 (*<n-l), x n = {-l) n+l xx + h i*» IT x 7- ( 42 - 10 ) 



X2>0,...,i/ n >0 2<p<n 

V-2 + ... + V n <d 



We will now fix the choice of the positive A in (|4.2.3p such that 



X<- -5- . (4.2.11) 

(n-l)d + r 

In the segment I° ut :— {e~ s < x n < 1} of W" (j — 1, 2), we choose qj — 1 different points Mji, . . . , M Jj9j ._i, 
and one point Mj qj 6 W" will be chosen in the segment — e~ s > x n > —1. Let Uj S denote the coordinate x n 
of Mj S (s = 1, . . . , qj). As NS = 1, near the segment I° ut the (N + l)-th iteration of the time-6 map Yg is the 
map T jA+5 from (l4~2~g|) . Thus, the map Y 5 Ar+1 near 7°"* will be given by 

5l = e ((n-2)A-l)i a .(n-l)A-l 42j = e -^ +1 (2<i<n-2), 

(4.2.12) 

Xn— i — Oj e x n , x n ^ ij e x n x\ 

(where we denote Xi = Xi at i ^ n— 1 and £E n _i = a; n -i ~bj). This map takes the segment J™* onto the segment 
{bj - 1 < x n - X < bj - e~ xs , xi = ... = x n - 2 = x n = 0} e W7. Let P j)S+1 = P,- w M js (s = 1, . . . , qj - 1), 
and let Pji be a point from {bj + er x& < x n _i < b 3 ■ + 1, x\ = . . . = x n _ 2 = x n = 0} € W?. By (|4.2.12j) . the 
coordinate x n -i of Pj, s +i equals to bj — e~ xs uj x . 

We take sufficiently large integer m and choose some points Pj s and Mj s , sufficiently close to Pj S and M/ s 
respectively (j = 1, 2; s = 1, . . . , qj), such that at s < qj — 1 we have Mj s = Lj m Pj s (where Lj t is the map 
(|4.2.7p ). At s = qj we assume Mj q . = Lj^n+^gP^. where lj is an integer to be defined later (see (|4.2.31[) ; note 
that IjS is uniformly bounded). Denote the coordinates of Pj s and Mj s as (z'j sl , . . . , z'j s n _ 2 , bj + z'j s n _ 1 , z'j sn ) 
and (u' sl , . . . , u' js „_ 2 ,bj + < 1; t/ sn ) respectively. By (|4.2.7p, 



- A ' m z; si (i = l,...,n-l), U ; sn = e m z; sn (4.2.13) 



at s < qj — 1. At s = qj we have 

^.i=e- Ai(m+ ^%^ (i=l,...,n-l), u' jq . n = e m+l > 5 z' qjn . (4.2.14) 

Note that u'j si are small at i < n — 1, as m is assumed to be large, and IjS is bounded. The values of z'- si with 
i =/= n — 1 will be taken sufficiently small as well, and we will keep 

u' jsn = u js and z' jiS+hn _ x = -e- x5 uj s x , (4.2.15) 

in order to ensure the closeness of Pj s to Pj S and Mj s to Mj S . 

The first of the small perturbations which we add to the map Yg is localized in a small neighborhood of 
the points Y s ~ (Pj jS +i) (so outside these small neighborhoods Yg remains unchanged). We take these localized 
perturbations such that in a sufficiently small neighborhood of Mj s = Y^ N (Y^ 1 Pj S+ i) the map Y g N+1 (where 
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Ys denotes the perturbed map) is given by (14.2. 12[) with the following correction term 

z > _ f_l)n+l.U-(n-2)A 



^n— 1 Oj 



Z js,n~l 



-Am / 



(4.2.16) 

+ ]^[ (ip — u jsp) P 

V2>0,'-,"n>0 2<p<Tl 
V2 + ... + V n <d 

added into the equation for x n , where Sj SL , are small coefficients to be determined later (see (|4.2.2ip ). The first 
term in (14.2. 16[) is small as well (see (|4.2.13l) . (14.2.14p 'l; in the second term the values of x n and z'j s n _ 1 are 
bounded away from zero, and the exponent ( j- — (n — 1)) is larger than r (see (|4.2.1ip ), hence the second term 
is also small with the derivatives up to the order r at least. The first two terms in (|4.2.16l) ensure, in particular, 
that at e = the coordinate x n of Y S N+1 M' JS coincides with that of Pj >s+1 (see (j4.2.13|) . (|4.2.15|) . (|IT2T3T) ). We 
want Pj. s+1 = Y s N+1 M' js at e = 0, so we put 



z j,s+l,l ~ e 11 js Z js2> 



(l-(n-2)A)«-Am,,(n-l)A-l / 

i'2 5 

(4.2.17) 

= e_A(Tn+i) 4.,i+i (2<*<^-2), 



(see (|4.2.12p . (|4.2.13p . (|4.2.15p ). At s = 1 we assume 



z' jU = at i<n-2, ^l.n-i = e" Xa/2 . (4.2.18) 

Now, the values of z^,i^ sl are defined by (|4.2. 13j> . (|4. 2. 14|) . f)4.2. 15]) . f|4.2. 17[) for all j, s,i. As one can see, z^ si 
at i ^ n — 1 and it^ si at i ^ n tend to zero as m —¥ +oo, i.e. Pj s — > P JS and Mj s — » Mj S indeed. 

At all small e the map f jtl+s oL jm = f d mN + N + 1 takes a small neighborhood of Pj s into a small neighborhood 
of Pj s+1 . We choose some T)(m) that tends to zero as m — ¥ +oo and some, independent of m, coefficients 
Cj S i > 0, and introduce rescaled coordinates v\, . . . , v n near Pj s by the rule 

xi\b d - x n -x\ n - l - l l x = 4 1 |^ S ,„_ 1 |- 1 - 1/A + C jsl r, e- Xm ^ Vl 

(4.2.19) 

4< = + C jsi n e-^-t-^Vi {2<i<n- 1), x n = z' jsn + C jsn 77 e~ m v n 

(recall that \bj — x n -i\ is close to 1 near Pj s , hence (I4.2.19|) is a smooth coordinate transformation). Since rj 
tends to zero as m —> +00, any bounded region of values of v corresponds to a small neighborhood of Pj 8 . 

After the rescaling, the map Tj,i+s ° Lj m = y™ n+n+1 from a small neighborhood of P!- s in to a small 
neighborhood of Pj s+1 takes the following form (see ^J^ , ^J^ , ^M , ^^M , ^IB ^ ^£^ ^ ^^M )' 

Cj tS+ i t iVi = e~ xs C J:Sti+ iv l+ i (i <n- 2), 

Cj <s +i, n -iv n -i = e~ XS (uj s x - (uj s + C jsn V v n )~ X )/r), 



Cj^ s + l, n V n — (— l) n+X (j)j s Cj s lV\ + ej su Ej sv Y\ 



i>2>Q,...,i<ri>0 2<p<n 

V-2 + ... + l> n <d 
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where we denote 

cf> ja = e s {u js +r ] C 3sn v n ) 2 -^ x \z l ^ n _ 1 +i 1 C J , s , n - 1 v n _ 1 \T- {n - 1 \ 

2<p<n 

Note that ^2<p<n-i( 71 — P) u p — {'n — 2)d, hence 1 — A Y^2<p<n-i( n ~ V) v v > (see (|4.2.11|1 ). Therefore, all the 
coefficients Ej SL , tend to infinity as m — > +00 (provided 77 tends to zero sufficiently slowly). 
As we see, by putting 



Cj,s+iA — e A<5 Cj, Sj i+i (i < n — 2), Cj>+i,n-i — ^ e XSu j S X ^^Cjsn, 

and 



(4.2.20) 



the map T^.i+a ° L jm near Pj s takes the form 

Vi = v l+1 (i<n-2), v n -i = v n + 0(77), 

5„ = (-If+V + 0(77) + £ fc iai/ [] ( 4 - 2 - 22 ) 

v 2 >o,...,^„>o 2<p<n 

i.e. it can be made as close as we want to the map Hj S , provided m is taken large enough (recall that 77 — > as 
?7i — » +00). We take 77 tending to zero sufficiently slowly, so, as we mentioned, £y„ — > 00 as m — > +00, which 
implies that all £j SU — > (see (|4.2.2f [) ), i.e. our perturbation to is arbitrarily small indeed. 

It follows that in the rescaled coordinates the map (Tj^i+s o Lj„^j = y^ 1 ( mjV +JV+i) f rom a small neigh- 
borhood of into a small neighborhood of P'^ q . can be made as close as we want to the map Hj q . o . . .0 Hji, 
provided m is large enough (the rescaled coordinates near Pjj and Pj q . are given by formulas (|4.2.19|) . where 
the coefficients Cju > are taken arbitrary, and the coefficients Cj qj i are recovered from the recursive formula 
(|4.2.20j) ; since m does not enter (|4.2.20[) . it follows that Cj q -i stay bounded away from zero and infinity as 

777 — > +0O). 

Now, from (|4.2.7p we obtain that the same holds true for the map Lj_ m+ i^o (Tj i + s o Lj„^j = Y^ 3+rnN+qj ^ mN+N+1 ^ 

from a small neighborhood of Pj x into a small neighborhood of M'j q . , where the rescaled coordinates (i>i, . . . , v n ) 
are introduced as follows: 

Xl \ Xn -i - = u'^lu^^X- 1 - 1 ^ + C ml 77 e-Kn-Dm-xi^ 

(4.2.23) 

Xi = U' jqji + C ml 77 e -*rn(n-i)-Xlj6 Vi ( 2 < i < n - 1), X„ = « iflJ + C Jw „ 7? e 1 ^, 

with the same constants C^w as above. 
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Recall that, by construction, the point Y^ +1 Mj lies in C/3-J.+ in the region a,s-j,+ + 1 — <5 < x n -\ < 
a3-j,+ + 1. We add to the map Y$ an additional perturbation, localized near the point YgMj q ., such that the 
corresponding map G^i+a = Y S N+1 will have the following form near M'j q ^: 

Xl — {Xnl 1 n 2 -* A (xi — e x ( m + l i s ) Zjgji (x n -\ — bj)/Zj q . n _ 1 ), 

i i = |a ! „| A (a!i-u^ < ) (2<*<n-2), ( 4 - 2 - 24 ) 
*n-i = 03-j,+ + 1 - ^ - In |z n |, x n = -|x„| A (x„_i - 6j- - u' jq .^_ x ). 

Note that Wj 9ji at i < n — 1 tend to zero as m — > +00 (see (|4.2.14[1 ). while the values of x n near and 
z j qj n-i are bounded away from zero; the exponent j — (n — 1) in the first line is larger than r (see (|4.2. 1 1[) ) . 
Thus, for sufficiently large m, map (|4.2.24j) is indeed a small perturbation of the map Gj^+s given by (|4.2.9|) . 

Denote P£_j = Gj t i+sM' jq .. By (|4.2.24|) . this is the point with the coordinates Xi = at i ^ n — 1, and 
x„_i = a 3 _j, + + 1 + (kj — 1)6 (we assume that the coordinate x n of Mj g . is u'j q . n — Uj qj = —e~ RjS where 
Kj E (0, 1] is defined by (|4.2.31[> ) . Introduce rescaled coordinates near P£_j by the rule 

x t = e-^+ K ^C jqji 77 e- Xm{ - n ^ Vl (i < n - 2), 

X n -1 = 03-j,+ + 1 + (1 - Kj)S + e^ + ^ s C ]q]n 77 V n - U (4.2.25) 
x = e~^ +K ^ sx C- i 77 e" Am ?j 

(with the same constants Cj g .j as above). In coordinates (I4.2.23|) . (I4.2.25P . map (|4. 2.241) takes the form 
(ui, . . . , u n _i,u n ) i — ^ . . . , u„, — V n _i) + 0(77), i.e. it becomes arbitrarily close to the map $0 (see (I3.23p ) as 
?7i — > +00. Thus, in the rescaled coordinates, the map Gj.i+goLj.m+i^o (T^i+iS o Lj m ^j = }>'j+( 1 + < ?j)( mA '+ Ar + 1 ) 
from a small neighborhood of P'^ x into a small neighborhood of P^-j q , can be made as close as we want to the 
map $0 o Hjg, o . . . o Hji as m grows. 

Analogously, we take the point Mj : {x„_i = dj_ — 1 + 6/2, Xi = (i ^ n — 1)} € Uj_, and perturb the 
map Ys near Y S N Pj in such a way that the map Qj,i+s = Y s +1 near Mj will be given by 

Xl = e ->>i(6-x»-i-i+*i-) Xi (i < n - 2), 

(4.2.26) 

r p-ACa-Xn-i-l+Oj-) - _ f-Xn-i-l+o^-- , p-m 

It is a small perturbation of the map Qj,i+s from (I4.2.6|) . and it takes Mj to P'^ (see (|4.2.13p . (|4.2.15p ). When 
we introduce rescaled variables near Mj by the rule 

Xi = e sx ^ 2 C jU 77 e-Amfn-i-l)^ (j < n _ 2), 

x„_! = 0j -_ - 1 + 5/2 + -e^C,-, i, n _i 77 «„_!, (4.2.27) 

^ ^ Gjin 77 6 7j n , 



map (|4.2.26[) will take the form v = v + 0(rj), i.e. it is close to the identity map. Thus, in the rescaled coordinates 



given by (14.2.271) . (I4. 2.25D . the map Gj,i+s ° Lj >m+ i jS o (fj,i+s ° L jm \ ' o Q j:1+S = Y t 



— ^lj+(l+q 3 )(rnN+N+l)+N+l 
r 8 
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from a small neighborhood of ML into a small neighborhood of P^_j , is as close as we want to the map 
<&o o Hj q . o . . . o Hji at m large enough, i.e. it is a close approximation of the map $j. 

Let us now determine the form of the map Sjt : v t— > v from a small neighborhood of Pj into a small 
neighborhood of Mj in the rescaled coordinates (I4.2.27l) . (l4.2.25p . By (14.2.51) . for an integer k > 0, the map 
Sj,ks = Y$ takes the point PL into ML if 



= (fe + /e 3 _.,- - 1/2)5-5 (4.2.28) 

(see (|4.2.27p . (|4.2.25p ). Since /3 -t +oo as fj, -)• +0 (see (|4.1.4p ). for every sufficiently large equation (|4. 2.281) 
has a solution jUj(fc), and Hj{k) — > +0 as /c — >• +oo. It follows that a(fij(k)) — > +oo. Thus, for any sufficiently 
large m we can find 7„± 6 (0, 1] and k such that 



— e ^3— j,g3-j ,n— 1 '/ e j 



(4.2.29) 



(4.2.30) 



e - 7n _«(w(fc)) = e - 5 / 2 Q ln 77 e" m . 

This guarantees that w„ = tpj(v n ) (see (|4.2.27p . (|4.2.25p . (|4T23|) V 

We also obtain Vi — vi at i < n — 2 by choosing 7y- g (0, 1] such that 

e - 7i+ «(^(fc)) = e a./2 Cjii 7? e -Am ; 

Finally, we fix the choice of the integer lj and kj G (0, 1] as follows: 

e^ +K ^ 5 = je x5 / 2 C 3 ^ Xn ^/C jqjn . (4.2.31) 

This (along with (|4.2.28[> ) gives us v n -i = v n -\ for the map Sj^s in the coordinates (|4.2.27|) . f|4.2.25|) . As wee 
see, the map Sj t ks i n the rescaled coordinates coincides with the map for v from some open neighborhood 
of D (if j = 1) or of $i o (if j = 2). 

Thus, we see that the map 

^2,1+4 ° L 2 .m+l 2 S ° (r 2 ,l+S ° L 2m ^ o Q 2 ,l+<5 ° S 2 .kS° 

° 61,1+4 o ^Tl,l+4 o Li m ^j o Ql,l+4 ° S\,kS = 

_ T>2fe+;i+i 2 + (2+gi+g 3 )(mJV+iV+l)+2(JV+l) 
— Z 4 

is a close approximation to the map F = $ 2 o $ 2 ° $1 o f 1, provided Hj q . o . . . o Hj\ are sufficiently close 
approximations to (j = 1, 2) and m is large enough. This completes the proof of the theorem. 



5 Birth of periodic spots from a heteroclinic cycle 

In this Section we finish the proof of Theorem [3] Let a C p -diffeomorphism / of a smooth two-dimensional 
manifold have a pair of saddle periodic points P and Q of periods p and, respectively, q. Denote as Tqi the map 
f p restricted onto a small neighborhood of P, and denote as Tq 2 the map f q restricted onto a small neighborhood 



2G 



n, + 



Figure 3: A non-transverse heteroclinic cycle 



of Q. By definition, XbiP = P and T^Q = Q- One can introduce coordinates (xi,y\) in the neighborhood of 
P and (^2,2/2) in the neighborhood of Q such that the maps T j will have the form 



where |Aj| < 1, \jj\ > 1; the dots stand for nonlinearities. The numbers Ai,7i and A 2 ,72 are the multipliers of 
the periodic points P and Q, respectively. We assume 



Every saddle periodic point lies in the intersection of two C p -smooth invariant manifolds: the points in the 
stable invariant manifold W s tend to the periodic orbit at forward iterations of the map, and the points in 
the unstable invariant manifold W u tend to the periodic orbit at backward iterations. With our choice of the 
coordinates, the unstable manifold is tangent at the periodic point to the y-axis, and the stable manifold is 
tangent to the a;-axis. One can locally straighten W u and W s , i.e. the coordinates (xj,yj) can be chosen in 



such a way that W^P) = { Vl = 0}, W^(P) = {x x = 0}, Wf oc (Q) = {y 2 = 0}, W%JQ) = {x 2 = 0}. 



Assume that W U {P) has an orbit of a transverse intersection with W S (Q). This means that := 
/ fel2 Mj~ G Wi oc (Q) for some positive integer fci2 and some point M 1 "(0,y 1 ") G W^ C (Q), and the curve 
/ fcl2 (Wj" c (P)) intersects W[ oc {Q) at the point M 2 H (a; 2 f ,0) transversely. We denote the map f kl2 restricted 
onto a small neighborhood of M-f as Ti 2 . It can be written as 



Xj — XjXj + . . . 



V] = IjVj + ■ ■ ■ . 



Ji:=|Ai7i|<l and J 2 := |A 2 7 2 | > 1. 



(5.1) 



X2 - x% = aixi +b 1 (y 1 -y 1 ) + 



y 2 = c\x\ +d 1 (yi-y 1 ) + ..., 



(5.2) 
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where e?i ^ because of the transversality of T 12 (W^ C (P)) to Wf oc (Q). As an orbit of intersection of W U (P) 
and W S (Q), the orbit Tpq of the point Mf is heteroclinic: it tends to the orbit of P at backward iterations 
of / and to the orbit of Q at forward iterations. Note that since d\ ^ 0, we may rewrite (|5.2[) in the so-called 
cross-form: 

d\{x 2 - x%) = Dxi + biy 2 + . ■ . , di{yx - y^) = y 2 - c\X\ + . . . , (5.3) 
where D := a\d\ — b\C\. 

Another assumption is that / has a heteroclinic orbit Tqp at the points of which W U (Q) has a tangency of 
order to with W S (P). This means that there exists a pair of points, M 2 ~(0,?/ 2 ~) € WjocW) and M{ l "(a;^,0) £ 
W)* (P), such that — f k21 M^ for some positive integer k 2 i, and the curve / (Wj^ c (Q)) has a tangency 
of order to with W Z * C (P) at (see Fig. 3). We denote the map f k21 restricted onto a small neighborhood of 
M% as T21. It can be written as 

xi - %i = a 2 x 2 + b 2 (y 2 - V2~) + ■ ■ ■ , Vi = C 2X 2 + d 2 (y 2 - y 2 ) m+1 + . . . , (5.4) 

where d 2 ^ 0; the dots stand for higher order terms. Obviously, to speak about the tangency of order to, the 
smoothness of / has to be sufficiently high, i.e. p > m + 1 . 

We will fix the orientation in the neighborhood of Q by requiring that the determinant D of the derivative 
matrix dTi 2 /d(x,y) is positive at the point M-j - . Then, we require that the determinant of dT 2 \/d(x, y) at 
the point is also positive, i.e. 62C2 < 0. This is always the case if the manifold on which / is defined is 
orientable and / is orientation-preserving. In Remark after Lemma [3] we discuss the case 62C2 > too. 

In Lemma [3] we will also need a technical assumption In I71I In I72I ^ In |Ai| In | A2 1 (one can always achieve 
this by an arbitrarily small perturbation of / without destroying the order m tangency between W U (Q) and 
W S {P)). In fact, we can always assume 

ln| 7 i|ln| 72 | <ln|A 1 |ln|A 2 |; (5.5) 

the case In | ^/i | In | -7-2 1 > In | Ai | In | A2 1 reduces to the given one by considering the map / _1 instead of / (and 
interchanging P with Q). 

Let us imbed / in any (m + l)-parameter family /„ of maps, of class C p with respect to coordinate and 
parameters. The corresponding map T%\ will also depend smoothly on parameters. Since it has form (|5.4|) at 
v = 0, i.e. all the derivatives of y~\ with respect to (y 2 — y 2 ) vanish up to the order to, it follows that at non-zero 
v the map T 2 \ can be written in the form 

m—l 

xi-xf = a 2 x 2 + b 2 (y 2 - y%) + . . . , f/i = c 2 x 2 + p, s (y 2 - y^Y + d 2 (y 2 - y 2 ) m+1 + ..., (5.6) 

s=0 

where [i s are smooth functions of v (note that x± and y^ also depend on v now: the value of y^ is fixed by 
the condition d m y~i/dy 2 n — at y 2 = y%\ thus y^ is C p-m -function of u, and fi s are also C p ~ m ; the high-order 
terms that are denoted by dots in (|5.6p and (|5.3p depend now also on v, C p_m -smoothly). 

Denote 8 — \ In Ji] In J±\. The value of 8 is a C p ~ 1 -smooth function of v. We may always put a family /„ in 
general position, i.e. we may further assume 

-^r, r- + 0. (5.7) 

This, in particular, means that by changing v we may change the values of any of the parameters 8 and \x s 
while keeping the other parameters constant. For example, one may change the value of 8 and keep /x = 0, i.e. 
keep the order to tangency between W U (Q) and W S (P). 
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The two periodic orbits (the orbit of P and the orbit of Q) and two heteroclinic orbits, Tpq and Tqp 
comprise a heteroclinic cycle. Let U be a small neighborhood of the heteroclinic cycle. Varying the parameters 
v can lead to the destruction of the heteroclinic cycle and to bifurcations of other orbits in U. We will further 
focus on one instance of such bifurcations. We call a periodic orbit of the map /„ single-round if it stays entirely 
in U and visits a small neighborhood of each of the points M^ 2 only once, i.e. the point of intersection of the 

single-round orbit with a small neighborhood of is a fixed point of the map T2iTq 2 2 Ti2X'o 1 1 for some positive 
integers k\ and k<i ; the image of this point by the map T*l lies in a small neighborhood of , the image by 
T\2TqI lies in a small neighborhood of M%, and the image by T^T^TqI lies in a small neighborhood of M% ■ 

Lemma 3 . There exist a sequence vi — > and sequences of integers ku — > +00 and k% — > +00 such that 
at v — vi the map /„ has a single-round periodic orbit which corresponds to a fixed point of the map 7/ := 
T 21 T^ l T 12 T^ 1 , and the map 71 near this point is, in some C p -coordinates (u,v), given by 

(u,v) =<t>( Ul v) + o(\u\ m + \v\ m ), (5.8) 

where $ denotes the time-1 map by the flow 

u = v, v = -$>(u)(l + v) (5.9) 

near (u,v) = (0,0); here ^ is a polynomial such that \l/(0) = and \&'(0) > 0. 

Remark. System l|5.9p has an integral H{u, v) = J ^(u)du + v - ln(l + v) = -£ + ^'(O) 2 ^ + . . .. Thus, when 
^'(0) > 0, the equilibrium at zero is a center: every orbit of (|5.9[) is in this case a closed curve surrounding 
the origin. The corresponding time-1 map $ is, therefore, conservative (it preserves the area form — du A dv) 
and has an elliptic point at the origin. Note that map (I5.8[) can, by an arbitrarily C" l -small perturbation, be 
made equal to (u,v) = $(u,v) identically in a sufficiently small neighborhood of zero. Hence, Lemma [3] im- 
plies, that by an arbitrarily C" l -small perturbation of the given map / a periodic point can be born in a small 
neighborhood of the heteroclinic cycle such that the first-return map near this point will be area-preserving 
and the point will be elliptic. We recall that one of the conditions of the lemma is that the heteroclinic cycle is 
"orientable" in the sense that the determinant of dT2i/d(x,y) at the point M 2 is positive. However, it is easy 
to show (see Lemma 8 in [S]) that a non-orientable heteroclinic cycle can always be perturbed in such a way 
that a new (double-round) orbit of heteroclinic tangency of order (to — 1) between W U (Q) and W S (P) is born, 
and the corresponding cycle is now orientable; applying the lemma to the newly born cycle, we find that an 
elliptic periodic point can be born by a perturbation which is arbitrarily small in the C m_1 -metric. In any case, 
we have that for any fixed r, given a C°°-map / which satisfies (|5 . 1 1) and which has a heteroclinic cycle with a 
sufficiently high order of tangency between W U (Q) and W S (P), an arbitrarily C r -small perturbation leads to 
the birth of an elliptic periodic point: a point whose both multipliers lie on the unit circle and, importantly, 
the corresponding first-return map is area-preserving. By |32j , an arbitrarily C""-small perturbation of the map 
in a neighborhood of such point creates an orbit of homoclinic tangency; the perturbation does not destroy the 
area-preserving property. By Theorem 5 of [6] , an arbitrarily C r -small perturbation of any area-preserving map 
with a homoclinic tangency leads to the birth of a periodic spot. In other words, once Lemma [3] provides us 
with a periodic orbit whose normal form is area-preserving up to a sufficiently high order, the birth of periodic 
spots from the heteroclinic cycles under consideration follows from the results of [B] on perturbations of area- 
preserving maps. 

Proof of Lemma [3 At to = 1 an equivalent statement can be found in [301 131] , so we further focus on the 
case to > 2. Since the periodic points P and Q are saddle, it follows that given any small Xjo and yjk (where 
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j = 1,2) and any k > there exist uniquely denned small Xjk and yjo such that (xjk,Ujk) = Tq (xjo, yjo) and 
all the points in the orbit {(£jo)> yjo),Toj(xjo,yjo), . . . , TJy (a^o , J/jo)} he in a small neighborhood of P (at j = 1) 
or Q (at j = 2), see [301 133] • Denote 

Xjk = \jXj + £jk(xjo, Vjk, v), yjo = lj k yjk + Vjkixjo, Vjk,v)- (5-10) 

By [33) , one can introduce C p -coordinates (xj , yj ) near the saddle periodic points in such a way that 

&* = o(Aj), %fe = o( 7 7 fc ), (5.11) 

i.e. the map TJy written in the "cross- form" (|5.10p is linear in the main order. By [7J, the same o(A^) and, 
resp., 0(77 ) estimates hold for the derivatives of the functions and r]jk up to the order (p — 1) with respect 
to (xjo,yjk) and up to the order (p — 2) with respect to parameters; while for the higher order derivatives up 
to the order p for which the number of differentiations with respect to the parameters does not exceed (p — 2) 
we have that they uniformly tend to zero as k — > oo. 

Let S > be sufficiently small. Consider rectangular neighborhoods LT^" : {\xj — Xj\ < S, \yj\ < §} and 
ELJ : {\yj — yj" | < 6, \xj\ < 5} of M- and MJ (j = 1,2), and take any sufficiently large integers ki and k 2 - By 

(|5.10p . (|5.1ip . the set = n^nT ~- 3 II7 is non-empty: it is a strip of the form {\"f kj yj — yj\ < 5 + o(l)/ s _ ) .+ 00 }. 
We introduce a new coordinate j/j on (Jjki such that 

2/j,oid = l~ k3 yj,new + iJ^Vjkj (Xj,yj, ne w) (5-12) 

(i.e. yj, n ew equals to y ]kj from (|5.10p ). 

By construction, the map T^T^ is defined on axki \ i n the new coordinates, when this map takes a point 
[xx,yx) G fiifej into a point (£2,2/2) G C2fc 2 it can be written in the following form (see (|5.3[) . (|5.12| 1): 



di(x 2 - 4) = DX^X! + fafa) + ba~ k2 y 2 + o(7 2 " fc2 ), 

di(yi - yi) = -c l \ kl x l - Mxi) + i2 k2 y2 + o{ l2 - k2 ), 



(5.13) 



where (f>x,2{xi) — o(X kl ). Analogously, the map T 2 iT 2 that takes a point (#2, 2/2) e cr 2fc 2 into a point (aci, j/i) G 
cr lfel is given by 

£i - = e^A^a^ + 62(2/2 - + °(#2 ~ V2) + (^2 2 )> 

m— 1 /_ < \ 

7i" fcl yi = caA* 3 ^ + £ - IfcT + - 2/2~) m+1 + - % + o^ 2 ) + o( 7l " fel ) " ' 

s=0 

(see (|5.6p ). We will now change the variable j/i to yj + d^ 1 (ciA5 :i a;i + 2 (^i))- Then the second lines in (|5. 13[) 
and (I5.14p will change to 

di(yi - vx) = i2 k2 m + o(j2 k2 ), 

and, respectively, 

m 

lx kl Vl = C2\ k 2 2 X2 + 5>.(lfi» - V2 Y + d 2 (y2 ~ V 2 ) m+1 + o{{ V 2 - V 2 ) m+1 ) + o{Xf) + ofrf*), 
s=0 

where ft s = fi s + 0(X kl ) (s = 0, . . . , m — 1) and p) m = 0(X kl ). 
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It is easy to see, that since d\ ^ 0, d ^ 0, one can find constants Cj = 0(\\3-j\ k3 - j + | 7 j| kj ) and 
Kj = 0(|73_j | _fc3 -j + |Aj| fcj ) (j — 1,2) such that after the following shift of the origin: 

Xj : new — Xj Xj &j ') Vj,new — ~l~ -^j'j (5.15) 



we will have (x 2 , n ew, Dx, new) = at (xi >n&w , y 2 , new ) = 0, andxi,^ = 0, d m yi, new /dy?; new = at (x 2 , new , V2,new 
0. Thus, after this transformation, the maps Ti 2 Tq^ and T 2 iTq 2 2 will be written as 

d lX2 = D\^ Xl + 0( l2 7 k2 y 2 ) + o((|A!| fcl + |7 2 r fc2 )a;i) 



k 2 „-. , „r„-k.-. \ , /. , I -ft, \ ( 5 - 16 ) 



and, respectively, 



dm = 7 2 ^2 + o( 72 fc2 y 2 ) + o(| 72 |- fc2 a:i), 
xi = b 2 y 2 + o{y 2 ) + 0(X k2 x 2 ), 



(5.17) 



7r fcl yi = c 2 \\ 2 x 2 + M + d 2 y? +1 + o(y? +1 ) + o(\ k2 x 2 ), 

s=0 

where the modified parameters jl s are such that fl s = /i s + o(l)/ Cl 2 ^+co- 

By virtue of (|5.7[) . by an arbitrarily small change of v we can make 8 :— |ln IA272I/ In |Ai7i| | rational without 
changing the values of fi. Therefore, we may from the very beginning assume that at v = the system has 
a heteroclinic cycle with a tangency of the order m, and the value 8q :— 8(0) is rational. Further we always 
assume that 

fci = 8 k 2 (5.18) 

in (I5.16[) . (|5.17l) . We will also assume that both k\ and k 2 are even, so \ kj and 7- 3 are positive. It follows from 
(15Tgl) . (l53j) .(l5Tl) that 

l 2 k2 < \ kl <C 7f fel < A2 2 • (5.19) 

Now, let us introduce new, rescaled coordinates (X\, Y±, X 2 , Y 2 ) and parameters (B,Eq, . . . ,_E m _i) by the 
following rule: 

7 j — fei/m -k 2 /m v -.fcj —kx/m -k 2 /m v 

x\=b 2 d 1 ^ l 7 2 Xi, z 2 = D&2A17! 7 2 A 2 , 

— ki/m — k 2 (l+l/m) v , —kx/m — k 2 /m v 

Vi=li 7 2 y i> y2=di7i 7 2 y 2, 

Ms=«i(7i 7 2 J ^ (s = 0,...,m-l), 

q = p + H-B/(Db 2 c 2 ) _ 
fe 2 |ln|Ai 7 i|| 

recall that Db 2 c 2 < by our assumptions, so the new parameter B should be positive. After the rescaling, 
the maps Ti 2 T l (given by (|5.16l0 and T 2 iT^ (given by (|5.17|l ) are rewritten as follows (we take into account 

X 2 = Xi + o(l)fc 1 2 _i. +00 

Yi =Y 2 + o(l) kl2 ^ +QC , and 

Xl =Y 2 + 0(l)fc! 2 ^ + 00, 

Y 1 = -BX 2 + Y^'o E s Yi + d + o(l) kl2 ^ +oc , 

where d := d 2 d™ +1 ^ 0. It follows that with our choice of k% <2 the map T := T 2i T 02 2 T i2 Tqi on ai kl can be 
written as 

m — 1 

X = Y + o(l) kl 2 ^ +00 , Y = -BX + E sY S + d Y m+1 + o(l) kl 2 ^ +oc , 

s=Q 
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where (X,Y) are the coordinates which were previously denoted as (Xi,Yi). By denoting the right-hand side 
of the equation for X as the new y-variable, we finally write the map T in the following form 

m— 1 

X = Y, Y = -BX+Y,E s Y s +dY m+1 +o(l) kl ^ +oc . (5.21) 

s=0 

The transformation to coordinates (X, Y) was of class C p with respect to the original coordinates and C p ~ m 
with respect to the parameters (more precisely: the m-th derivatives with respect to the coordinates are C p ~ m 
with respect to both the coordinates and parameters). The o(l)-terms in (|5.2ip are functions of (X,Y) and 
(£?, Eq, . . . , E m —±) which tend to zero as fci.2 — > +00 along such derivatives up to the order p, for which the 
number of differentiations with respect to the parameters does not exceed (p — m). In any case, since p > m + 1, 
we have that their derivatives with respect to (X, Y) up to the order m tend to zero, all with at least one 
derivative with respect to (B, E , . . . , E m -i). 

We mention that a similar form for the rescaled first-return map was obtained in [3T] for the case m = 1. 
Map (|5.21[) was also obtained in [B] as the rescaled first-return map near a homoclinic tangency of order m. 
Let us show that if the o(l)-terms are sufficiently small, then a map of form (15.211) has, at certain uniformly 
bounded values of E and B > and in a bounded region of (X, Y), a fixed point near which the map is smoothly 
conjugate to (|5.8p . This will give us the lemma: by (|5.20l) . (|5.15[) . any bounded values of (X,Y) correspond to 
(xi,yi) — > (xi ,Ui) as fei ( 2 — > +00, i.e. if the values of (X,Y) are bounded at the fixed point, then the latter 
belongs to the domain of definition of the first-return map ^iT^r^T^ 1 at all sufficiently large fci,2, and any 
bounded values of (Eq, . . . , E m -i) and B > correspond to v k — > as fei^ — > +00. 

We may take any Y* as the coordinate X = Y = Y* of the fixed point of (|5.2ip : then the second equation 
defines the corresponding value of E = Y*(l + B) - J^L\ E s(Y*) s ~ d (Y*) m+1 + o(l). After shifting the 
origin to the fixed point, the map takes the form 

Y = -BX + 2Y + JXi E S Y S + 0(y m +!) + 0(\XY\ + X 2 ) ■ (l)^ +oo , (5 ' 22) 
where the coefficients E a , s = 1, . . . , m — 1, are related to E m —i by a diffcomorphism: 

m—l 

&i = Ei-2+Y, jE^Yy- 1 + (m + l)d (Y*) m + o(l), 



3=2 



m — l 



and 



E S = E S + E j C s j (Yy- s + C s m+1 d(Y*r+ l - s +o(l)(s = 2,...,m-l), 

3=s+l 

E m = (m + l)dY* + o(l), B = B + o{l). 



As we see, bounded values of B, E\, . . . , E m give us bounded values of B, Eq, . . . , B OT _i as well. 

We will further fix B = 1 and will take (E±, . . . , E m ) = o(l)fc^ +oc , so the map will limit, as k — > +00, to 

X = Y, Y = -X + 2Y. 
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The fixed point of this map at (X, Y) = has a double multiplier (+1). Let us recall the normal form theory 
for maps with such a fixed point. In general, we can always write such map in the form 



u = u + v+ ^2 a lj u l v j + o(\u\' m + \v\ m ), 



2<i+j<m 



2<i+j<m 



P ijU V +o(\u\ m + \v\ m ). 



(5.23) 



By denoting z — v + ^2<i+j<m atijuV + o(|u| m + \v\ m ) = u — u, the map takes the form 



z+ Y ^uV+o(H m + \z\ m ). (5.24) 



2<i+j<m 



Here 



<f>ij = fcj + J2 C k a k,i+i-k + Gij, (5.25) 

i<k<i-\-j 

where Gij stand for polynomial functions of the coefficients on*ji and fiyji with i' + j' < i + j- each monomial 
in has degree 2 or higher. After the polynomial coordinate transformation (u, z) n- (U, Z), where 

U = u+ Aju n - j z\ Z = z+ A j (u n - j z : > -u n - j z j ), 

0<j<n-2 0<j<n—2 

the map will retain its form (|5.24j) with the coefficients 4>ij unchanged for i + j < n, and 4>ij^ new = 4>ij + 
So<J<i— 2^ ~ l ~ ^)^n-!^ at i + j = n. As we see, by choosing Ai appropriately, one can make 4>ij with 
j > 2, i = n — j assume any given values. Thus, by making such transformations consecutively, with n running 
from 2 to to, one can make all the coefficients 0y with j > 2 vanish; the new coefficients <f>ij (j — 0, 1) differ from 
the old ones by nonlinear (i.e. quadratic and higher order) expressions which depend only on the coefficients 
4>i>ji with i' + j' < i + j , i.e. they are still given by (15.251) with some new Gij. Obviously, for any fixed to, we 
arrive to a similar conclusion for any map (with a fixed point at zero) whose linear part is sufficiently close to 
that of (|5.23p . Specifically, for any sufficiently small £1.2, a map of the form 



u = u + v+ y a i jU % v i + o(\u\ m + \v\ m ), 

2<i+j<m 

v = (l + e 1 )v + e 2 u+ J2 +o(\u\ m + \v\ m ), 

2<i+j<m 



(5.26) 



can be brought by a smooth coordinate transformation to the normal form 

m 

u = u + z, z=(l + e 1 )z + e 2 u + y((j) i QU i + (j) l - hl u''z) + o{\u\ m + \z\ m ), (5.27) 



i=2 

where 



<fao = Pio + Gia + HiQ, <t>n = pix + {i + l)a i+Xfi + G iX + H iX . (5.28) 

Here is a polynomial of a^ji and fli'ji with %' + j' < i + j, all monomials of Gij are of degree at least 2; the 
function Hij is a linear function of ayji and ft/j/ with i' + j' = i+ j; the coefficients of Gij and Hij depend on 
£1,2 as well, and Hij vanish at e\ = £2 = 0. 
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In fact, given any m, for any sufficiently small ex,%i there exists a pair of functions -0o,i( M ) sucn that after 
an appropriate choice of the coordinates (u, v), map (|5.26l) can be made o(|u| m + |v| m )-close to the time-1 map 
by a flow of the form 

u = v, v = ipo (u) + vtpi (u) . (5.29) 
Obviously, the flow has to have an equilibrium at zero and the eigenvalues of the corresponding linearization 

„ ^ |, the linearization matrix of (15.261) . Thus, 

£2 l + £l / ' 



matrix have to be logarithms of the eigenvalues of 
we assume 



V>o(0)=0, V'i(O) =l n(l+ £l -g 2 ), 
^(0) = - ln(I + | + \Js2 + ^) HI + | - 



£2 



F 2 

-)• 

4 ; 



(5.30) 



The time-i map of (|5.29[) can be found by expanding it in powers of the initial conditions (u(0),v(Q)). This 
gives us the time-1 map in the form which is brought, by an 0(ei j 2)-close to identity linear transformation, to 
form ([5~2f)|) with 



+3,0 




C 



i+3 



•0(ei 



+ 0(£i, 2 ) 



+ tpi+j-1,1 

c, 



c 



i+j-l 

j 



0(e x 



3 



+ 0(£i, 2 ) 



(5.31) 



where onj and /?y are nonlinear functions of the coefficients ?/Vj' with i'+j' < i+j (we denote ipj(u) = tpijU 1 ). 
Two maps of form (|5.26[) can be made o(|w| m + |?;| m )-close by means of a smooth coordinate transformation if 
their normal forms ()5.27|) coincide up to the order m. Thus, we find from formulas (|5.28[) . (|5.30p . (|5.31|) that a 
map of form (|5.26|) can indeed, after an appropriate coordinate transformation, be made o(|w| m 4- |i>| m )-close to 
the time-1 map by a flow of form (|5.29[) with the coefficients defined by 



tpio = PiO + i>- 



(0, 



4>ii = fti + (i + l)(a i+ i ) o - ft+i,o) + i>u, 



(5.32) 



where ipy is a polynomial function of a^j' and /3i>ji with i' + j' < i+j; the monomials that include a^y or fy/ji 
with i' + j' < i + j are of degree 2 or higher, while the terms that include a^j' or fy/j/ with i' + j' = i + j are 
linear and the corresponding coefficients are small of order 0{e\^)- It follows from (|5.30j) that formulas (|5.32|) 
remain valid at j — 1, if we formally put /?io = £2, An = £1, c^io — a oi = 0. 

Returning to map (|5.22p . by putting X = u and Y = u + v, we bring the map at B = 1 to form 

u = u + v, v = v + p (u) + vPi(u) + 0(v 2 ) + o{\u\ m + \v\ m ), 

where the coefficients of the polynomials (3j = Y^i=i J Pij u% U = 0, 1) satisfy 

/?i0 = Ei + o(l) fc _>. +00 , fin = (i + l)E i+ i + o(l) fc ^+oc- 

e sufficiently small, make the 



Thus, by making normal form transformations described above, we can, at E\ 
map o(|u|™ 1 + |i>| m )-close to the time-1 map of the flow (|5 . 29[) with 



iO 



Ei + o(l)k^+oo + o(Ei, . . . , E m ), Vil = o(l)fc^+oo + o(Ei 



■ E m )', 



see (|5.32p . It follows that at sufficiently large k we can always choose the values of E\ , ... , E m in such a way 
that tp (u) = sipi(u), where s = 1 if ^(0) < 0, and s = -1 if ip[(0) > 0. 
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As we see, at an appropriate choice of the parameters the map can be made o(|u| m + |i;| m )-close to the 
time-1 map of the flow 

u = v, v = ipi(u)(s + v), 

which takes the desired form (15.91) after the change (u, v) —t (su, sv); here ty(u) — —ipi(su), so the requirement 
^'(O) > is ensured by our choice of s. □ 

This work was supported by grant ISF 273/07 and by Royal Society grant "Homoclinic bifurcations", and it was 
done during my stay at MPIM, Bonn. 
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